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A MATHEMATICAL PROCEDURE FOR PREDICTING THE TOUCHDOWN 

DYNAMICS OF A SOFT-LANDING VEHICLE 

By George A. Zupp, Jr., and Harold H. Doiron 
Manned Spacecraft Center 

SUMMARY 


A mathematical procedure for predicting the touchdown dynamics of a soft-landing 
vehicle is presented. The vehicle is idealized as a pseudorigid structure to which a 
landing gear system is attached. The landing gear system is made up of telescoping 
struts which house shock-absorbing mechanisms and connect the vehicle structure with 
footpads that contact the landing surface. The forces acting on the idealized vehicle 
during landing are from the vehicle landing gear system, the reaction control system, 
the descent-stage rocket engine thrust, rocket-engine-nozzle-crushing loads, and grav- 
ity. Forces acting on footpads are due to forces from the landing gear struts and from 
the landing surface. The elastic properties of the vehicle structure and of the landing 
gear system are approximated by assuming that the structure deforms linearly with the 
applied load. Correlation between analytical and experimental data is presented for a 
l/6-scale lunar module model. Conclusions from these data are presented. 


INTRODUCTION 


A primary objective of the Apollo Program is to soft-land a manned spacecraft, 
the lunar module (LM), on the lunar surface. During the first critical seconds of touch- 
down, the landing gear system must absorb the kinetic and potential energies of the 
spacecraft without causing the LM to topple and must attenuate the landing loads to pre- 
clude damage to the spacecraft structure or possible injury to the astronauts. To en- 
sure that the landing gear system has the capability to arrest the spacecraft motion 
safely, landing gear performance must be established. Although methods of determin- 
ing landing gear performance for aircraft systems are well established, these methods 
are not directly applicable to spacecraft, primarily because of the differences in con- 
figurations. Consequently, development of new analytical and experimental procedures 
has been necessary to determine LM landing performance, as well as to aid in the de- 
velopment of the landing gear system. One of the more elaborate analytical procedures 
that has been developed for simulation of touchdown dynamics and for landing gear de- 
velopment is discussed in this report. Although developed for the LM, the procedure 
is applicable to other soft-landing spacecraft. 

In the initial development phase of the Apollo LM landing gear system, analytical 
simulations of touchdown dynamics by a NASA contractor (ref. 1) and by the NASA 



Manned Spacecraft Center (MSC) (ref. 2) were primarily restricted to landings on rigid 
surfaces that resulted in planar motion. These t 3 ^es of landings were considered crit- 
ical for vehicle stability and landing gear energy absorption. Analysis of the touchdown 
dynamics of rigid models at the NASA Langley Research Center (refs. 3 and 4) indi- 
cated that, from a stability standpoint, nonplanar landings could be more severe than 
planar landings. As a result of this finding, more emphasis was placed on the simula- 
tion of nonplanar landings for LM touchdown analysis; consequently, the touchdown 
analysis presented in this paper evolved. 

The purpose of this report is twofold; 

1. To document the analytical procedures used in developing a mathematical 
model in the LM lunar landing dynamics 

2. To make these procedures available to others who may be confronted with the 
task of mathematically modeling the touchdown dynamics of a soft-landing vehicle 

This report includes a definition of pertinent coordinate systems and coordinate trans- 
formations. The equations of motion governing both rigid- vehicle and landing-gear- 
footpad dynamics are presented, together with the numerical procedure employed to 
solve the equations. Considerable detail is provided in defining the kinematic relation- 
ships used to model the articulating landing gear geometry. Various forces generated 
in the landing gear struts and the manner in which these forces are computed are dis- 
cussed. A general description of the techniques used to model a reaction control sys- 
tem (RCS), descent engine rocket thrust, and descent engine nozzle crushing is also 
provided. Both rigid and nonrigid landing surface models are discussed, and consid- 
erable detail is provided on the manner in which various surface irregularities and 
footpad constraints are simulated. Data from the digital simulation analysis and a 
1/6-scale LM model drop-test program are compared, and conclusions regarding the 
accuracy of the simulation technique are made. 

The authors are particularly indebted to Joseph A. Reid, formerly of ITT Federal 
Electric Corp. , for his efforts in the assembly and preparation of this report. The 
authors also wish to express sincere appreciation to John A. Schliesing and Ben W. 
Holder, MSC, for their technical comments and suggestions. 


SYMBOLS 


AMP amplification factor for the descent propulsion system 

rocket thrust due to the nozzle exit plane approaching 
the landing surface 


AP 


area of the bottom of a footpad, in 


2 


APF 


frontal area of a footpad in contact with soil, in 


2 


2 



A . 

n,j 


AB 


ACM 


A 

AI . 

n,] 


th 

instantaneous vector directed from the n " footpad center to 

th 

a fixed point on the j surface curb 

vector directed from footpad A to footpad B 

vector connecting the rigid-vehicle center of mass and a 
landing gear footpad 

unit vector normal to the j " surface curb and directed 

th 

from the initial surface contact point of the n " footpad 


a . = DWALLI . cos WALANG . 
n, j n, j n, J 


B1 

n 


B2 ,B3 ,B4 
n’ n’ n 


vector directed from the origin of the body coordinate sys- 
tem to the primary-to-secondary-strut attachment point 
th 

of the n ^ landing gear assembly 

vectors directed from the origin of the body coordinate sys- 
tem to the rigid -structure attachment points of the pri- 

th 

mary and secondary struts in the n^^ landing gear 
assembly 


Ki 


th th 

unit vector directed along the j surface curb of the n"^” 

footpad 


BP 


n 


vector directed from the origin of the body coordinate sys- 

th 

tern to the footpad of the n “ landing gear assembly 


b . = DWALLI . sin WALANG . 
n,J n,] n,] 


CLUEW 

n 

CP 


cs 


ce ,ce ,C9 

x’ y’ z 

Cco , Cw , Cw 
x’ y’ z 


computer program input clue used to define the curbs on the 

th 

landing surface for the n^^ footpad 

constant used to compute primary strut stroking-velocity- 

dependent force, (lb/(ft/sec)) ^ 

constant used to compute secondary strut stroking -velocity- 
's 

dependent force, (lb/(ft/sec)) 

commanded values for the Euler angles used in the RCS 
simulation, deg 

commanded attitude rates used in the RCS simulation, 
deg/ sec 


3 



D2 

n 


3 

constant used in numerical integration (eq. (20)), ft/sec 
til 

length of the primary strut inner cylinder, ft 

til 

initial bearing separation in the n*^ primary strut, ft 


DIST 


DWALL . 
n,3 


DWALLI . 
n, j 


desired maxim.um allowable movement per integration step 
between the footpads and the rigid-vehicle center of 
mass, ft 

th 

instantaneous distance between the n footpad center and 
ttl 

the j " constraining curb, ft 

normal distance between the n^^ footpad center and the 
curb when the footpad initially contacts the surface, ft 


D0 ,D6 ,D0 
x’ y’ 


z 


El 


deviations in Euler angles from the commanded values, deg 

2 

constant bending stiffness of a primary strut, Ib-ft 


EX, EY, EZ 


combined attitude and attitude rate errors used in RCS 
simulation 


EXM, EYM, EZM 

E0 ,E0 ,E0 
X y z 

Eo) , Eu) , Ew 
X y z 


^B, x’ ^B,y’ ^B,z 


deadband widths for RCS simulation 

attitude errors for the body axes used in RCS simulation, 
deg 

attitude rate errors used in RCS simulation, deg/sec 
magnitude of descent engine thrust at touchdown, lb 
components of sloshing force vector (eq. (115)), lb 




force per unit vertical tank acceleration caused by a unit 
sloshing amplitude, (lb/ft)/(ft/sec^) 

total force applied to the tank wall by the sloshing propel- 
lant, lb 

force at bearing 1 of the primary strut (fig. 14), lb 
force at bearing 2 of the primary strut (fig. 14), lb 
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F31 ,F41 
n n 

FBR 

FBRn 

FC31 ,FC41 
n n 


FCP^ 


instantaneous axial force in the secondary struts of the n*"” 
landing gear assembly, lb 

friction forces present during strut stroking, lb 

th 

primary strut bearing friction force in the n*"" primary 
strut, lb 

secondary strut stroke-dependent forces in the n'"" landing 
gear assembly, lb 

th 

stroke-dependent force in the n " primary strut, lb 


FPl, FP2 

FRCSZ 
FSCl, FSC2 

FSTl, FST2 

FV31 ,FV41 
n’ n 

FVP 

FVP^ 

FX, FY, FZ 

FXP ,FYP ,FZP 
n n n 

FZI 

n 


first- and second-stage crushing forces of the honeycomb 
shock absorbers in the primary struts, lb 


Z -component of FRCS, lb 

compression stroking forces of the secondary struts for 
the first- and second-stage honeycomb cartridges, lb 

tension stroking forces of the secondary struts for the 
first- and second-stage honeycomb cartridges, lb 

th 

stroking-velocity -dependent forces in the n^^ landing gear 
secondary struts, lb 

viscous forces present during strut stroking, lb 

th 

stroking-velocity -dependent force in the n*^” landing gear 
primary strut, lb 

forces on the idealized rigid vehicle resolved along the X-, 
Y-, and Z-axes, respectively, lb 

th 

forces on the n " footpad in the X-, Y-, and Z-directions, 
which are independent of the instantaneous sliding direc- 
tion of the footpad, lb 

force component normal to the landing surface applied to 
th 

the n " footpad by the landing gear, lb 
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F 

F* 

TF 

F^ 

n 


FCMG 

FDRAG 

n 

FDRAGP 

n 


FFL 


FGRND 

n 


FM 

n 

n 

FP 

n 


FPI 


FRCS 


FRCSI 

FS 

n 

FT FI 


total force vector acting on the idealized -vehicle center of 
mass 

total force vector due to engine thrust and the nozzle failing 
load 

instantaneous axial force in the primary strut inner cylinder 
th 

of the n*^^ landing gear assembly, produced by compres- 
sion stroking of the primary strut 

total landing gear forces acting on the vehicle center of mass 
expressed in the inertial coordinate system 

th 

soil drag force acting on the n ^ footpad 

th 

frictional drag force on the bottom of the n " footpad 


nozzle failure load vector 
th 

drag force on the n^^ footpad with direction opposite to the 
instantaneous sliding direction of the footpad in the 
X-Y plane 


th 

force normal to the primary strut at the n footpad 

component of the secondary strut forces normal to the 

th . . . 

n primary strut 

resultant force on the n*'^ footpad due to landing gear forces 

resultant of the n^^ landing gear forces acting on the rigid- 
vehicle center of mass, expressed in the inertial coordi- 
nate system 

net force vector produced by the RCS jets 

RCS jet force in the inertial coordinate system 

vector sum of the secondary strut forces acting on the 
n primary strut 

force vector on the rigid-vehicle center of mass due to en- 
gine thrust and nozzle failing load, expressed in the in- 
ertial coordinate system 
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FTHRUST 


FWF„ , 
n, j 


thrust force of the descent propulsion system rocket 

friction force between the n*^ footpad and the surface 
curb 


FXYP 

n 


th 

force vector on the n " footpad independent of the instantan- 
eous sliding direction of the footpad 


GOXX, -GOXY, GOXZ, components of the body coordinate system to inertial coor- 

GOYX, -GOYY, -GOYZ, dinate system transformation matrix (eq. (2)) 

-GOZX, GOZY, GOZZ 

G gravitational acceleration vector expressed in inertial 

coordinates 


gjj’ Sy’ ^2 

HOLE 

n 

I ,I ,I 
XX yy zz 

^xx, f’ ^yy, f’ ^zz, f 
I 

xy 

I 

xz 

I 

yz 

ij.k.n 

K 


components of the gravitational acceleration vector ex- 
pressed in the inertial coordinate system along the X-, 

2 

Y-, and Z-axes, respectively, ft/sec 

parameter which defines the surface elevation for the 
n^^ footpad, ft 

mass moments of inertia about the X„-, Y„-, and Z„-axes, 

2 D D a 

respectively, slug-ft 

modified mass moments of inertia in the body coordinate 
system of the rigid vehicle, corrected for exclusion of 

2 

the sloshing mass, slug-ft 

cross product of inertia with respect to the X„-Y„ plane, 

9 a D 

slug-ft^ 

cross product of inertia with respect to the X_-Z_ plane, 

9 ±5X5 

Slug-ft 

cross product of inertia with respect to the Y„-Z„ plane, 

9 ±5 ±5 

slug-ft^ 

indices 

constant characterizing the thrust tailoff properties, 1/sec 
bearing strength of the soil at the surface, psi 




increase in soil bearing strength with depth, psi/ft 
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^3 

Increase in soil bearing strength with footpad velocity, 
psi/(ft/sec) 

^4 

increase in soil bearing strength with the square of the foot- 


2 

pad velocity, psi/(ft/sec) 

^5 

static soil resistance to horizontal movement, psi 

^6 

increase in soil resistance to horizontal motion with re- 
spect to footpad lateral displacement, psi/ft 

^7 

increase in lateral resistance of the soil with respect to 
the depth of penetration of the footpad, psi/ft 

^8 

increase in soil resistance with respect to the square of 

2 

the footpad horizontal velocity, psi/(ft/sec) 


spring constant representative of the elastic property of 
the honeycomb shock-absorber material, Ib/ft 

^P, B 

primary strut backoff spring rate without trunnion-point 
elasticity, Ib/ft 


primary strut onset spring rate without trunnion-point 
elasticity, Ib/ft 


spring constant representative of the slope of the load- 
deflection curve of the spacecraft structure at a trunnion 
point, Ib/ft 

^SC, B 

secondary strut backoff spring rate in compression stroking 
without trunnion-point elasticity, Ib/ft 

^SC,o 

secondary strut onset spring rate in compression stroking 
without trunnion-point elasticity, Ib/ft 

^ST, B 

secondary strut backoff spring rate in tension stroking with- 
out trunnion -point elasticity, Ib/ft 

^ST, 0 

secondary strut onset spring rate in tension stroking without 
trunnion-point elasticity, Ib/ft 

K,pp 

trunnion -point spring rate at the primary strut- structure 
attachment point, Ib/ft 


trunnion-point spring rate at the secondary strut-structure 
attachment point, Ib/ft 
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*^P,B 




^SC,B 


^SC,o 


^ST, B 


^ST,o 


KB3^, KB4^ 


Ic k" Ic k Ic Ic 

11’ 12’ 21’ 22’ 31’ 32 


modified backoff compression spring rate for the primary 
struts, Ib/ft 

modified onset compression spring rate for the primary 
struts, Ib/ft 

modified backoff compression spring rate for the secondary 
struts, Ib/ft 

modified onset compression spring rate for the secondary 
struts, Ib/ft 

modified backoff tension spring rate for the secondary 
struts, Ib/ft 

modified onset tension spring rate for the secondary struts, 
Ib/ft 

spring rates for secondary struts as a result of primary 
strut bending, Ib/ft 

constants which model a given RCS control mode 


L21 


n 


L31 ,L41 ,L2P 
n’ n’ n 


M 

^ 0 ’ ^ 1 ’ ^2 
Mj 

M„ 


th 

distance along the n " primary strut from the upper attach- 
ment point to the secondary strut attachment point, ft 

undeformed lengths of the secondary and primary struts of 
th 

the n landing gear assembly, ft 

any given distance below the quiescent propellant surface 
(fig. 28), ft 

distance from the spacecraft center of mass to the plane of 
the quiescent propellant surface (fig. 29), ft 

distance from the spacecraft center of mass with propellant 
to the spacecraft center of mass without propellant 
(fig. 29), ft 

mass of the idealized rigid vehicle, slugs 

constants characterizing the rocket engine thrust tailoff 
properties 

fixed mass of residue propellant, slugs 

sloshing mass of residue propellant, slugs 
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M(s) bending moment as a function of length along the primary 

strut, ft-lb 

^ 

MRCS RCS torque vector 


MS 

n 


th 

moment vector perpendicular to the n primary strut, 
caused by action of secondary struts 


MTF 


moment vector acting on the idealized vehicle due to the 
rocket thrust and nozzle failure load 


m 

OBST 

n 


PFRC . 
n, J 

PMASS 

PR 


number of landing gears 

parameter which defines the height of constraining obstacles 
for the n*^^ footpad, ft 

til 

point where the primary strut is joined to the n " footpad 

th th 

crushing force of the n footpad against the j surface 
curb, lb 

effective mass of each footpad, slugs 
undeformed footpad radius, ft 


PSTRK , 
n,J 




Rl, R2,R3 


R, 


Rn 


th 

variable indicating the amount of n ^ footpad crushing normal 
to the curb, ft 

sloshing amplitude, ft 

rigid -vehicle acceleration of the propellant tank normal to 

2 

the vertical axis of the tank, ft/sec 

distances from the center line of the LM defining the landing 
gear geometry (fig. 23), ft 

vector from the origin of the body coordinate system to the 
centroid of the nozzle failure load 

vector from the origin of the body coordinate system to the 
centroid of the nozzle exit plane 


S2P 

n 

S31 ,S41 
n’ n 


th 

compression stroke of the n primary strut, ft 

th 

stroke of the secondary struts in the n ** landing gear assem- 
bly, ft 
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lateral translation of a footpad, ft 

vector directed from the rigid-vehicle center of mass 

normal to a vertical plane through two adjacent footpads 

position coordinate along the length of the primary strut, ft 

transformation used to transform body angular rates to 
Euler angular rates 

transformation used to transform vector components from 
body system coordinates to inertial system coordinates 

transformation used to transform Euler angular rates to 
body angular rates 

transformation used to transform vector components from 
gravity system coordinates to inertial system coordi- 
nates 

transformation used to transform vector components in the 
inertial system coordinates to the body system coordi- 
nates 

transformation used to transform vector components from 
inertial system coordinates to gravity system coordi- 
nates 


*^f x’ ”^1 y torques about the Xg- and Yg-axes, respectively, due to 

propellant sloshing, ft- lb 

T , T ,T resultant torques about the X_-, Y„-, and Z„- axes, 

^ y ^ respectively, ft-lb « B B 


TORQX 


magnitude of the torque generated about the X„-axis by 
RCS thrusters, ft-lb ^ 


TRCSX, TRCSY, TRCSZ 
T 

TAF 


TCMG 


RCS torques about the X„- , Y_-, and Z„-axes, respec- 
tively, ft-lb a a a 

net torque about the idealized-vehicle center of mass 

total engine thrust force, the product of the amplification 



factor AMP and the normal thrust force FTHRUST 

resultant torque about the rigid-vehicle center of mass 
due to landing gear strut loads 


t 


time, sec 


11 


p 

u. 


n 


UW 


n 


V, 


V, 


V 

VEL 


V 


V21„ 


V2P^,V31_^,V41„ 



WALANG 


time of LM descent engine thrust termination command, sec 

time at which descent engine thrust starts to tailoff, sec 

til 

time at the beginning of the integration step, sec 

reference start time for engine thrust tailoff curve, sec 

time at which the LM footpad touches the surface, sec 

til 

coefficient of sliding friction between the n footpad and 
the landing surface 

til 

coefficient of friction between the n*"" footpad and surface 
curbs 

horizontal touchdown velocity, ft/sec 


overturning velocity, which is a measure of the degree of 
instability of an unstable landing, ft/sec 

vertical touchdown velocity, ft/sec 


absolute value of the maximum relative velocity between 
the rigid-vehicle center of mass and the landing gear 
footpads, ft/sec 

translation velocity vector of the rigid-vehicle center of 
mass 


th 

vector directed along the center line of the n^^ primary 
stnit from the upper end of the strut to the secondary 
strut attachment point 

til 

vectors defining the instantaneous strut geometry of the n^” 
landing gear assembly (fig. 5(b)) 

position vector locating the intersection point of 

DWALLI . with a curb 
n,J 

angle from the positive X-axis to a line in the X-Y plane 
normal to a curb, deg 


X, Y,Z 


inertial coordinate system axes 
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body coordinate system axes 






gravity coordinate system axes 


XB4 , YB4 , ZB4 ,XB3 , X_, Y„, and Z_ coordinates of several points defining 

n n n n d Jd 

YB3^, ZB3^, XB2^, YB2^, landing gear geometry (eq. (26)) 

ZB2 ,XB1 YBl ,ZB1 , 
n n, n n 

XBPj^, YBPj^, ZBP^ 

th 

XP , YP , ZP instantaneous position coordinates of the n^ footpad in the 

" “ inertial coordinate system, ft 


XPI ,YPI ,ZPI 
n n n 

Zl, Z2, Z3, ZC 
ZNOZ 


th 

inertial coordinates of the n*^ footpad when it first contacts 
the landing surface, ft 

spacecraft geometrical parameters, ft 

surface clearance of the centroid of the rocket-engine- 
nozzle exit plane, ft 


ZPEN' 

n 


vertical penetration of the n 


th 


landing gear footpad, ft 


a 


slope of the landing surface, deg 




angle defining secondary strut geometry (fig. 23(b)), deg 


y. 


s 


exponents of the primary and secondary strut stroking ve- 
locities used to compute stroking-velocity-dependent 
forces 


At 


numerical integration step size, sec 


^3, n’ ®4, n 


bending deflections of the primary strut due to secondary 
strut loading, ft 


deflection of a secondary strut trunnion point due to com- 
pression loads (fig. 11(b)), ft 


6 . 

i,n 


apparent change in length of the i “ secondary stimt due to 
th 

n primary strut bending, ft 

deflection of the primary strut trunnion point (fig. 11(a)), ft 


deflection of a secondary strut trunnion point due to tension 
loads (fig. 11(b)), ft 
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n 

bending deflection of the center line of an idealized primary 
strut, ft 


th 

bending deflection of the n “ primary strut at the secondary- 

n 

to-primary- strut attachment point due to strut loading, 
ft 

0 ,0 ,0 
x’ y’ z 

Euler angles, deg 


primary strut bearing friction coefficient 

“n 

natural sloshing frequency, rad/sec 

0)„} 6D 

angular rate vector components of the rigid vehicle, ex- 
pressed in the body coordinate system, deg/sec 

X y z 

2 

(a 

sloshing frequency squared per unit vertical tank acceler- 
ation (eq. (109)), (rad/sec^)/(ft/sec^) 

Operators: 

(•) 

first derivative with respect to time 

C‘) 

second derivative with respect to time 

C) 

vector quantity 

n 

first derivative of a vector with respect to time 

C) 

second derivative of a vector with respect to time 

A 

( ) 

unit vector 

[ ] 

matrix 

[ ]-' 

matrix inverse 

[ 1^ 

matrix or column vector transpose 


Symbols that represent vectors are also used without the vector notation in this 
paper; omission of the vector notation indicates absolute value of the vector. 
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GENERAL DISCUSSION OFTHE MATHEMATICAL MODEL 


The formulation of the mathematical model for predicting touchdown dynamics 
of a soft-landing system starts with the definition of the governing dynamics equations. 
In this case, the soft-lander cab is given three degrees of freedom in translation and 
three degrees of freedom in rotation. Each landing gear footpad is allowed three de- 
grees of freedom in translation. The governing dynamics equations are the well-known 
Euler's equations for rotation and Newton's equations for translation. A numerical 
integration technique is employed to obtain the solution of these dynamics equations. 

The significant forces acting on the idealized soft-landing vehicle are assumed 
to be those from the landing gear struts, the thrusting rocket engine, the RCS, the 
collapsing rocket engine nozzle, and gravity. The forces acting on the idealized foot- 
pads are due to surface-footpad interaction and forces from the landing gear struts. 
Propellant sloshing forces are also usually present during touchdown. When the pro- 
pellant sloshing masses are small, the sloshing forces are considered to be negligible, 
and the propellant mass is treated as a rigid part of the vehicle structure. (In the 
appendix, a linear sloshing model is developed. ) Care must be exercised in the appli- 
cation of linear sloshing models to touchdown dynamics problems because of the large 
sloshing amplitudes experienced during some touchdown conditions. 

The landing gear is idealized as an arrangement of tubular struts with telescop- 
ing members, and the strut forces are assumed to be a function of Coulomb friction, 
strut stroke j and strut stroking rate. The landing gear and vehicle structural elasticity 
is approximated by the analogy that the structure deforms linearly with an applied load. 
However, the present model has been adapted for use with a more complex dynamic 
elastic model of the spacecraft structure which is described in reference 5. 

The RCS and engine thrusting forces are guidance and control parameters and 
are functions of the landing procedure. A general discussion of the methods by which 
these forces are modeled is included. 

Forces on the footpads result from landing gear strut forces and from resisting 
forces generated by the landing surface. Various landing surface models are dis- 
cussed. In the early stages of the Apollo Program, the uncertainty of the composition 
of the lunar surface, from a bearing capacity standpoint, resulted in considerable work 
by a NASA contractor to define realistic soil surface models. The soil models devel- 
oped by this contractor were incorporated into the mathematical touchdown simulation 
described in this report and were used to evaluate footpad penetrations for LM landings 
on postulated lunar surfaces (ref. 6). The soil model discussed in this report is a 
general model that allows the prediction of soil-footpad forces as a function of footpad 
movement and velocity. A rigid landing surface is also considered where the footpad 
is constrained upon surface contact or is allowed to slide under the influence of fric- 
tion and constraining curbs. 
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COORDINATE SYSTEMS AND TRANSFORMATIONS 


Three basic coordinate systems are used in the formulation of the mathematical 
simulation of the touchdown dynamics. Each coordinate system is an orthogonal right- 
hand system. 


Inertial Coordinate System (X-, Y-, and Z-Axes) 

The inertial coordinate system is defined such that the X-Y plane forms the land- 
ii^ surface. The Z-axis is perpendicular to and directed up from the X-Y plane. All 
of the rigid-vehicle and landing-gear-footpad equations of motion are expressed in the 
inertial coordinate system. 


Gravity Coordinate System (X^-, Y^-, and Z^-Axes) 


In the gravity coordinate system, the negative Z^-axis coincides with the gravity 
vector, and the X^-axis coincides with the inertial system X-axis. The inertial sys- 
tem and the gravity system are related as shown in figure 1. The angle a between 

the Y^- and Y-axes is the slope of the 


Gravity vector 

I 

I 



landing surface. The coordinate trans- 
formation [Tj^] from the inertial sys- 
tem to the gravity system is given by 


Pio] = 


10 O' 

0 cos a sin a 
0 -sin a cos a 


( 1 ) 


Figure 1. - Gravity and inertial co- 
ordinate system orientation. 
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Since equation (1) is an orthogonal transformation matrix, the transformation fTp,] 
from the gravity system to the inertial system is *■ -* 

'l 0 0 ~ 

[■^Gl] = f Ig] = [^Ig] ® 

_0 sin a cos a_ 

Initial-touchdown-velocity input data are expressed in the gravity system and then 
transformed to the inertial system by the transformation [TqjI • 


Body Coordinate System (Xg-, Yg-, and Zg-Axes) 

The body coordinate system is fixed in the landing vehicle so that the origin coin- 
cides with the idealized- rigid- vehicle center of mass, and the Zg-axis is directed up- 
ward parallel to the vehicle vertical center line (fig. 2). The body system is, related 
to the inertial system by the set of Euler angles 9 , d , and 9 illustrated in fig- 

3 X V z 



(a) Top view. (b) Side view. 

Figure 2. - Orientation of the body coordinate system in the rigid vehicle. 
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(a) Body and inertial systems coinciding. 


(b) 


Rotation about X_-axis 


(0 ). 



(c) Rotation about Y^’-axis (0^). (d) Rotation about Zg”-axis (0^) 


Z 



(e) Final orientation. 
Figure 3. - Euler angle definition. 



The rotation order is 0 , 9 , and 9 , with each successive rotation taken about 

X y z 

the indicated body axis in its position after the previous rotation. For example, the 
two coordinate systems can be oriented so that they coincide as shown in figure 3(a). 
The first rotation 9 is about the X„-axis as shown in figure 3(b). The next rotation 

X o 

9 is about Y„', the new Y_-axis position, as illustrated in figure 3(c). For the 

y o 

final rotation, the angle 9 is rotated about the new Z„-axis position, as in- 

Z x3 J3 

dicated in figure 3(d). The final orientation of the body coordinate system with respect 
to the inertial coordinate system is shown in figure 3(e). These rotations can be con- 
sidered to be a pitch- roll-yaw sequence for a pilot facing in the direction of the posi- 
tive Yg-axis. 

The transformation [TgJ which transforms vector components in the body 
system to vector components in the inertial system is 







GOXX 

-GOXY 

GOXZ 





II 

GOYX 

-GOYY 

-GOYZ 






_-GOZX 

GOZY 

GOZZ_ 

where GOXX 


COS 

9 

cos 9 








y 

z 





GOXY 

— 

cos 

9 

sin 9 








y 

z 





GOXZ 

_ 

sin 

9 









y 






GOYX 

- 

sin 

9 

sin 9 cos 9 

+ cos 9 

sin 

9 





X 

y 

Z X 


z 


GOYY 

- 

sin 

9 

sin 9 sin 9 

- cos 9 

cos 

9 





X 

y z X 


z 


GOYZ 


sin 

9 

cos 9 








X 

y 





GOZX 


cos 

9 

sin 9 cos 9 

- sin 9 

sin 

9 





X 

y 

Z X 


z 


GOZY 


cos 

9 

sin 9 sin 9 

+ sin 9 

cos 

9 





X 

y 

Z X 

z 


GOZZ 


cos 

9 

cos 9 








X 

V 






( 3 ) 


This transformation is also orthogonal, and jYjgJ which transforms inertial system 
coordinates to body system coordinates is given by 




(4) 
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Another important matrix transformation, which is employed, relates the rigid- 
vehicle angular rate vector components expressed in the body system to Euler angular 

rate components. The Euler angular rate components (0 , 0 , and 0 ) are directed 

X y z 

along the axis about which the respective Euler rotation (0 , 0 , or 0 ) is made. This 

X y z 

transformation is not orthogonal because the Euler angular rate components are not 
mutually orthogonal. 

The transformation from body angular rates to Euler angular rates is given by 


pBE_ 


cos 0 
2 

COS 0 


sin 0 


sin 0 
z 

cos 0 


cos 0 


sin 0 cos 0 sin 0 sin 0 
y z y z 


cos 0 


cos 0 


y ■ 


( 5 ) 


The transformation from Euler angular rates to body angular rates is given by 


_"^EB_ ['^BE_ 


-1 


Thus 


cos 0 cos 0 

sin 0 

0 

y z 

z 


-cos 0 sin 0 

cos 0 

0 

y z 

z 


sin 0 

0 

1 


( 6 ) 


0 



X 


X 

0 

y 

= Pbe] 

03 

y 

0 


CO 

z_ 


_ z 


( 7 ) 


where 0^, 0^, 0^ = Euler angular rate components 

0 ) , 0 ) , a> = body angular rate components 
X y z 
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A time derivative of equation (7) yields an expression for the time rate of change 
of Euler angular rate components. 


9 


db 


03 

X 


X 


X 

9 

y 

Pbe] 

“y 

w) 

03 

y 

9 

z_ 


S. 


“z. 


( 8 ) 


For programing purposes, equation (8) can be reduced to the following more efficient 
form: 


05 COS 9 - (h sin 6 - 29 6 + 0 o) 
e = X zy z zyyz 

X cos 9 


9 = db sin 9 + Jb cos 9 +99 cos 9 
yx zy zxz y 


9 = -9 sin 9 -99 cos 9 + 
z X y X y y z 


(9a) 

(9b) 

(9c) 


Equations (9a) to (9c) express the angular acceleration components cb , 

X 


db , and 


CO 


in a form that can be numerically integrated to yield the angular coordinates 9 , 

X 

and 9^ (the three Euler angles), which are used to describe the orientation of the 
rigid vehicle in the inertial coordinate system. 

A more common procedure for computing Euler angles is to integrate cb 


and CD to obtain co , 

Z X 


CD , and CD and then to transform these angular rates by 

y ^ 




means of equation (7) to obtain equations for Euler angular rates (for example, as in 
reference 4). The Euler angular rates are then integrated to obtain Euler angles. 
However, this approach is not used in the mathematical procedure developed in this 
report. The formulation given by equations (9a) to (9c) is more convenient for the par- 
ticular numerical integration procedure used to solve the equations of motion. This 
numerical integration procedure is discussed in a subsequent section of this report. 


EQUATIONS OF MOTION 


The landing vehicle is idealized as a rigid body of constant mass to which landing 
gear assemblies that have identical geometry and mass characteristics are attached. 
The mass of each landing gear assembly is assumed to be a point mass concentrated 
at the footpad of the assembly. Instead of using the actual mass of the landing gear 
assembly as the footpad mass, an "effective footpad mass" is computed so that the 
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effective footpad mass located at the footpad position approximates the inertia charac- 
teristics of the landing assembly. The effect of this approximation on the simulated 
landing dynamics has been studied in depth in reference 7. The conclusions in refer- 
ence 7 imply that if the mass ratio of the vehicle landing gear assembly to the vehicle 
is small (approximately 1 : 50 for the LM), the effect of this approximation on the sim- 
ulated landing performance is negligible. 

The instantaneous position coordinates of the rigid- vehicle center of mass (X, 

Y, and Z) together with the instantaneous values of the Euler angles {d , 9 , and 6 ) 

X y z 

specifying the orientation of the rigid vehicle, and the instantaneous position coordi- 
th 

nates of the n^" landing gear footpad XP , YP , and ZP in the inertial coordinate 

system allow for determination of the instantaneous geometry of the n*'“ landing gear 
assembly. As will be explained in more detail in the following sections, a knowledge 
of instantaneous landing gear geometry allows the computation of the instantaneous 
forces acting on the rigid vehicle and on each of the landing gear footpads. Therefore, 
the equations of motion which must be solved in the simulation of the landing impact 
dynamics are the three translational equations of motion of the rigid-vehicle center of 
mass, the three rotational equations of motion of the idealized rigid vehicle, and the 
three translational equations of motion for each landing gear footpad. 

The three translational equations of motion of the rigid-vehicle center of mass 
are obtained by the summation of all forces acting on the rigid vehicle and by the appli- 
cation of Newton's Laws of Motion. 


X 


Y 


Z 


FX 

M ^x 


FY 

M 


FZ 

M ^z 


( 10 ) 

( 11 ) 

(12) 


where FX, FY, FZ = summation of forces on the idealized rigid vehicle resolved along 

the X-, Y-, and Z-axes, respectively 

M = mass of the idealized rigid vehicle 

X, Y, Z = inertial accelerations of rigid-vehicle center of mass 


^x’ ^y’ ^z 


components of the gravitational acceleration vector expressed in 
the inertial coordinate system 


The three rotational equations of motion express the time rate of change of the 
angular momentum vector of the rigid spacecraft in terms of Euler angles and Euler 
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angular rates. These equations of motion given as follows were developed in the pre- 
vious section. 


(h . COS 6_ - d).. sin 6_ - 26_9_ + 6_ 




e = 

X 


z y z z y y z 
cos 9 


(9a) 


9 = cIj sin 9 + CO cos 9 +99 cos 9 
y X z y z z X y 


9 = -9 sin 9 -99 cos 9 + co 
z X y X y y z 


(9b) 

(9c) 


where iy, and are the components of the angular acceleration vector ex- 

pressed in the body coordinate system. These components are determined by solving 
the following Euler moment equations for rigid bodies (ref. 8). 


T=I 0)-I d)-I db+ (-1 CD -I CD +I CD^CD 

X XX X xy y xz z \ xy x yz y zz z) y 

- f-1 CD + I CD - I CD ICD 

\ xy X yy y yz zj z 

T =-I CD +I CO CD +(I CD “I CD “I CD JCD 

y xy^x yy y yz z V xx x xy y yz zy z 

+ (I CD +I CD "I CD ICO 

V xz X yz z zz zj y 

T=-I CD -I CD+I CD+ (~1 CD + I CD “I CD ^CD 

z xz X yz y zz z V xy x yy y yz z) x 

+ (-1 CD + I CD + I CD "\CD 

V XX X xy y xz z) y 


(13) 


(14) 


(15) 


where cu , cu , cd 
X y z 

I ,I ,I 
xx^ yy’ zz 

I 

xy 

I 

xz 


angular rate vector components of the rigid vehicle expressed in 
the body coordinate system 

mass moments of inertia about the and Z^-axes, 

respectively 

cross product of inertia with respect to the X„-Y„ plane 

o r> 

cross product of inertia with respect to the X^-Z„ plane 

13 13 


23 



= cross product of inertia with respect to the Yg-Z^ plane 
Ty, = sum of the torques about the Xg-, Yg-, and Zg-axes, respectively 

The footpad equations of motion describe translational footpad motion only and 
are obtained by summing all forces acting on the effective footpad masses and by apply- 
ing Newton' s Laws of Motion. The equations have the general form given in the follow- 
ing three equations, although the actual equations solved may vary with the specific type 
of landing surface postulated. The actual equations used for various types of landing 
surfaces are explained in detail in later sections. 


. I FGRND XP 

^Pn = PM^ F^n- . „ " ^ 1-^gx 




. I FGRND YP 

YP„ = /FYP_ - " - 1 + g„ 

v2 y 


n PMASS n 


y(xP„f,(YPn)= 


ZP = 


n PMASS 


1 (fzp \ 

ASS nj 


+ g„ 


(16) 


(17) 


( 18 ) 


th 

where XP , YP , ZP = instantaneous position coordinates of the n*"” footpad in the 

^ inertial coordinate system 

PMASS = effective mass of each footpad 

th 

FGRND = magnitude of a drag force on the n^ footpad with direction 

opposite to the instantaneous sliding direction of the footpad in 
the X-Y plane 

FXP^, FYP , FZP = forces on the n footpad in the X-, Y-, and Z-directions, re- 
spectively, which are independent of the instantaneous sliding 
direction of the footpad 

gx> g j gg = components of the gravitational acceleration vector expressed 
^ y ^ in the inertial coordinate system 


The equations of motion (eqs. (9a) to (12) and (16) to (18)) are solved by a numer- 
ical integration procedure. The initial conditions for these equations are specified in 
the following manner. The spacecraft attitude at the instant of impact is specified by 
the three Euler armies 0^, 0^, and 9^. These angles determine which footpad is ini- 
tially in contact with the landing surface, and this initial point of contact is taken to be 
the origin of the inertial coordinate system (0, 0, 0). The position of each footpad in the 
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body coordinate system prior to impact can be specified by known vehicle geometry. 
Let the vector expressed in body coordinates and directed from the origin of the body 

coordinate system to the footpad in contact with the surface be denoted by BP^. Then, 

the vector joining the origin of the inertial coordinate system and the origin of the body 
coordinate system is given by the vector transformation 


X 

Y 

Z 



(19) 


since the tip of is located at inertial system coordinates (0, 0, 0). The three ve- 

locity components of the rigid-vehicle center of mass are specified in the gravity coor- 
dinate system and are then transformed into inertial system components by means of 
the transformation given in equation (2). Thus, the three initial inertial coordinates of 
the rigid-vehicle center of mass given by equation (19), together with the three velocity 
components at impact, comprise the initial conditions for the numerical integration of 
equations (10) to (12). 

The angular rate vector components of the spacecraft at impact are initially spec- 
ified in the body coordinate system and are then transformed to Euler angular rates by 
equation (7). These Euler rates, together with the initially specified Euler angles, pro- 
vide initial conditions for the integration of equations (9a) to (9c). Equations (13) to (15) 
are solved as a set of simultaneous algebraic equations for the quantities cby, 

and db at each integration step prior to the numerical integration of equations (9a) 
z 

to (9c). After equations (9a) to (9c) are integrated, new values of co , , and 

X y z 

are computed by transforming the new Euler rates by means of the transformation 
given in equation (6). 

The positions and velocities of the landing gear footpads are determined by inte- 

th. 

gration of equations (16) to (18) only if the n'^ footpad is in contact with the surface. If 
th 

the n*"" footpad is not in contact with the surface, the position and velocity of the footpad 
are computed by assuming the footpad to be rigidly connected to the vehicle. The initial 
conditions for equations (16) to (18) are taken to be the inertial position and velocity of 
each footpad at the instant the Z- coordinate of the respective footpad is computed to be 
less than or equal to zero. If a footpad breaks contact with the surface, the landing 
gear geometry is assumed to remain unchanged from the instant the footpad leaves the 
surface until it recontacts the surface, and during the period of no contact with the sur- 
face, the position and velocity of the footpad are computed as though the footpad were 
rigidly connected to the vehicle. At the instant of recontact with the surface, the foot- 
pad equations of motion are initialized as before. 
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NUMERICAL INTEGRATION OFTHE EQUATIONS OF MOTION 


It has been mentioned that the equations of motion are solved by a numerical inte- 
gration procedure. There are many factors which should be considered in the selection 
of a particular method of numerical integration. Among these factors are the accuracy, 
speed, stability, integration- step- size control, and ease of adaptation of the method to 
the particular equations to be integrated. Among the more sophisticated methods now 
popular are (1) predictor- corrector methods that use backward difference formulas and 
that maintain a specified accuracy in the integration and (2) the well-known Runge-Kutta 
formulas. The predictor- corrector methods are cumbersome to use in problems where 
frequent changes in the integration step size are desirable. In attempting to maintain 
a specified accuracy, these methods will often perform an integration step forward, in- 
cur excessive integration error, and then waste the associated computations by making 
the same integration forward with a smaller step size. It was reasoned that because of 
the large number of computations involving the determination of landing gear geometry 
and forces necessary for each integration step, wasted integration steps could lead to 
unjustifiable amounts of computer time for landing simulations. Consequently, the 
predictor- corrector method was not used. To use a fourth-order Runge-Kutta formula 
for the integration of the rigid-vehicle equations of motion, it would be necessary to 
evaluate the landing gear geometry and associated force determination equations four 
times for each integration step. Since this large number of evaluations would also lead 
to very large amounts of computer time for the landing simulations, Runge-Kutta for- 
mulas were not used. A compromise integration scheme which was finally selected 
sacrifices a guaranteed numerical accuracy for a considerable gain in integration 
speed, while maintaining the desirable feature of variable step- size control. 

The selection of the numerical integration method presented in the following para- 
graphs was influenced by several factors which are peculiar to the landing simulation. 
First, it was recognized that small integration steps would be necessary because of the 
abrupt on-off nature of some of the forces which occur during landing impact, such as 
those obtained from RCS jets, multiple landing gear impacts, and rocket engine nozzle 
impact with the landing surface. Second, if integration were performed with relatively 
large integration steps and if a footpad impacted the surface during an integration step, 
then the integration step would have to be discarded and attempted again with a smaller 
step size. Another characteristic of the forces experienced by the vehicle during land- 
ing is that except for the discrete occurrences of abrupt force level changes, there are 
intermediate periods where the force levels remain fairly constant. 

These characteristics suggested that a simple method which assumed constant 
forces (or acceleration) over short integration steps could be expected to give reason- 
ably accurate results. To accommodate those integration steps where abrupt force 
chaises take place, it was assumed instead that the acceleration on the center of mass 
is a linear function of time during the integration step. 

Taking equation (10) as an example, a linear variation of acceleration over time 
step At implies 





( 20 ) 
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th 

where tj^ is the time at the beginning of the k integration step, and c is a constant 

for the integration step. If At = t - t , then the acceleration during the time step is 
described by ^ 


X = X. + c(t - tj^) 




A first integration of equation (21) from tj^ to tj^^j = tj^ + At yields 


Xf^Af = X. +X. At +i c(At)^ 

tj^+At tj^ tj^ 2 


( 21 ) 


( 22 ) 


A second integration of equation (21) yields 


\+At = \ * + 5 * I 

An estimate of X and X at t^ + At is obtained from equations (22) and (23) by 

assuming c = 0. With these estimates, new vehicle geometry, forces, and accelera- 
tions are computed, and the coefficient c is evaluated by computing FX/M at tj^ + At 
and then using equation (20) to obtain 


c 


~ M^Af 


(24) 


Then, equations (22) and (23) are used to obtain corrected values for X and X. A 
similar scheme is used for the integration of equations (9a) to (9c), (11), and (12). The 
footpad equations of motion are integrated by a standard fourth-order Runge-Kutta 
method where the magnitude of all forces on the footpads is assumed to be constant 
over an integration step, but where the variation in direction of the sliding velocity of 
the footpad during the integration step is taken into account. 

The logic used to control integration step size for this particular touchdown dy- 
namics problem is based on an ’’engineering decision, ” rather than on the usual trun- 
cation error checks used for determining step size for many variable step integration 
methods. This engineering decision is founded on the belief that as long as the amount 
of landing gear stroke or footpad movement is kept small (approximately 0. 1 inch for 
an LM landing) during an integration step, the forces on the vehicle and the footpads 
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will not vary svifficiently so as to violate significantly the assumption of linear varia- 
tion in acceleration over an integration time step. If the desired maximum allowable 
movement between footpads and the rigid-vehicle center of mass per integration step 
is specified by the input value DIST, then at any point in the integration, the integra- 
tion step size for the next step At is computed from the relation 


At = 


DIST 

VEL 


(25) 


where VEL is the absolute value of the maximum relative velocity between the rigid- 
vehicle center of mass and the landing gear footpads. The value for At computed by 
equation (25) is used, unless this value exceeds a specified maximum integration step 
size or is less than a specified minimum step size. 


IDEALIZATION OF LANDING GEAR GEOMETRY 


It has been mentioned previously that once the instantaneous forces acting on the 
rigid vehicle are determined, solution of the equations of motion is possible. The most 
significant forces which act on the rigid spacecraft during touchdown are generated in 
the landing gear struts. Therefore, detailed representation of these strut forces is 
necessary to produce an accurate simulation of the vehicle landing dynamics. To rep- 
resent landing gear strut forces accurately, a detailed simulation of landing gear artic- 
ulation during impact is necessary. 


The mathematical model of the soft-landing gear system considers two basic 
types of landing gear assemblies. One type is the cantilever gear presently used for 
the LM landing gear system (fig. 4(a)), and the other is the inverted tripod gear 


(fig. 4(b)). (The term ’’cantilever” in the 



(a) Cantilever (b) Inverted tripod 

assembly. assembly. 


Figure 4. - Two types of landing gear 
assemblies. 


name ’’cantilever gear” is incorrect in that 
the gear is not a cantilever in the geomet- 
ric sense. The designers of this landing 
gear originally referred to it as the ’’bend- 
ing gear” because of the bending of the 
primary strut under secondary strut load- 
ing. The name ’’bending gear” was not 
selected because it erroneously implied 
that strut bending was a primary function 
of the landing gear system. ) 

The landing gear strut assembly in 
either of the two arrangements is composed 
of one primary strut and two secondary 
struts. Each secondary and primary strut 
is composed of an outer cylinder and an 
inner cylinder, with the inner cylinder 
telescoping into the outer cylinder. In the 
cantilever arrangement, the lower end of 
the secondary struts is attached to the 
lower end of the primary strut outer 
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cylinder. For the tripod assembly, the lower ends of the secondary struts are attached 
to the lower end of the primary strut inner cylinder. Since the mathematical simulation 
of the tripodal arrangement is a special case of the cantilever arrangement, the follow- 
ing analysis of the landing gear geometry is developed for the more general cantilever 
configuration. A typical cantilever assembly is depicted in figure 5 by a set of vectors. 



(a) Orientation of landing gear assem- 
bly in body coordinate system. 


The landing gear represented in 
figure 5 is applicable to each of n land- 
ing gear assemblies; hence, n = 1, 2, 

3, .... All numbered points in figure 5(a) 
are assumed to be ball-joint connections. 
Points 2, 3, and 4 are trunnion points 
where the primary strut and the secondary 
struts are joined to the idealized rigid 
structure. In the present discussion, 
these points are considered to be fixed 
points on the body. The point P^, where 

the primary strut is joined to the footpad, 

can move in any direction as a result of 

stroking of primary and secondary struts. 

If the coordinates of P are known con- 

n 

tinually during a touchdown simulation, 
then vector equations can be constructed 
to define the instantaneous strut geometry, 
and subsequently, strut stroking velocities 
and strut lengths can be determined. 

Every point in figure 5(a) can be 
represented by a position vector joining 
that point with the origin of the body coor- 
dinate system. Then the vector equiva- 
lents of points 4, 2, 3, 1, and P^^ are 

B4 , B2 , ^ , fit , and ^ , respec- 
n n n n n 

lively, with body coordinates defined by 



Figure 5. - Landing gear strut vector 
diagram. 


B^n = [XB4„, YB4„, ZB4jT ' 

“n = [XB2„,YB2^,ZB2jT 

5S„ = [XB3„,YB3__,ZB3j'^ 
Mn = [XBl„.YBl„,ZBqT 


(26) 
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For the present, assume that the components of BP^ are known. (The solution of the 

footpad equations of motion will yield the position vector BP^. ) The vectors B2^, 

B3 , and B4 are known from geometric definition of the vehicle. The vector between 
n’ n 5 ^ ° 

points 2 and P is V2P , where 
n n 

(YBP„ - YB2„), ^BP„ - ZB2 jl (27) 


The modulus of the vector V2P , denoted by V2P , determines the instantaneous 

n’ n’ 

length of the primary strut. The vector V2P^ unitized is given by 


V2P 


^ 


n 


n 


(28) 


The vector 



is determined from the vector equation 


B1 = V21 + B2 (29) 

n n n 

where 

V 2 I = V^P ('L21 \ (30) 

n n \ n/ ' ' 


and the distance L21 between points 2 and 1 is constant. In the tripodal arrangement, 

^ — \ — X 

point 1 coincides with point P^; therefore, Bl^ is equal to BP^, and equations (29) 
and (30) are not necessary. (Of course, in the tripodal arrangement, L21^ is not 
constant. ) The vectors joining points 3 and 4 to point 1 are computed from 


V31 

= B1 

- B3 

n 

n 

I 



^ 



V41 

= B1 

- B4 


n n 


n 


(31) 
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The lengths of vectors V31^ and V41^ yield the instantaneous secondary strut lengths. 
The following equations define the unit vectors of equations (31). 


V31 


V41 


V«„'l 

n - V31_, I 

n - V41_^ J 


(32) 


The amount of stroke in each strut of the landing gear is required for determining 
axial forces in the struts. Figure 6 shows the primary and secondary struts in stroked 

and unstroked positions. If the undeformed 
lengths of the primary and secondary 
struts are denoted by L2P and by L31 

and L41^, respectively, the instantaneous 

strokes of each are given by 


Undeformed 
primary strut 


Compression 

stroke 


■ L2Pn- 


ft-. • 

V?Pn 

r 

vtrn 




KU S2Pf;. 


n-| 




S2P 


V2P - L2P 


n 


(33) 


(a) Primary strut. 


and 


-V3U 


Tension stroke 


Undeformed secondary stroke 


-L31n 


[Ib 


Compression stroke 


1 


H k K 
-i \-S31n 


S31 = V31 - L31 

n n n 


S41 = V41 - L41 
n n n, 


(34) 


(b) Secondary strut. 

Figure 6. - Stroked and unstroked pri- 
mary and secondary struts. 


The compression stroke of a secondary 
strut is a negative quantity, and a tension 
stroke is a positive quantity. 


FORCES IN THE LANDING GEAR STRUTS 


The problem of accounting for axial forces in the landing gear struts is rather 
involved because of force contributions from various sources. The various sources 
contributing to the landing gear strut forces are analyzed separately, summed, and 
resolved at the landing gear footpad. 
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Crushable-Type Shock Absorbers 

The primary type of shock-attenuation mechanism for the LM landing gear system 
is a sequence of crushable metal honeycomb cartridges housed inside the landing gear 
struts. These cartridges collapse at a constant load when the landing gear struts are 
stroked. Several honeycomb strengths may be stacked in series to form a desired 
load -stroke output. The present analysis provides for a maximum of three cartridges 
in the primary struts and four cartridges in each secondary strut. The primary strut 
has shock-attenuation capability in compression stroking only, and the secondary struts 
have shock-attenuation capability in both tension and compression stroking. Figure 7 
illustrates a typical honeycomb configuration for a primary and secondary strut 
assembly. 



Figure 7. - Primary and secondary strut 
crushable shock absorbers. 


Honeycomb material can be designed 
to fail at a constant load as the honeycomb 
material is crushed. If this constant load 
is predetermined, characteristic relation- 
ships between strut axial loads and strut 
strokes can be described by a load-stroke 
diagram. A typical load-stroke diagram 
for the primary strut is shown in figure 8. 



Figure 8. - Primary strut load- stroke 
diagram. 


The primary struts have only compression- type crushing of the honeycomb car- 
tridges. Tension loads are due to friction. The honeycomb cartridge loads are mod- 
eled in the following manner. In figure 8, the unstroked primary strut begins at 
station 0 and undergoes compression stroking. The strut load will increase linearly 
with stroke to the station A. The plateau denoted between stations A and B represents 
the failure load of the weakest or first-stage honeycomb cartridge. When the primary 
strut is stroked to station B, the first stage of the honeycomb cartridge has been ex- 
hausted. As the primary strut is stroked past station B, the load will again build up 
linearly with stroke vintil the second plateau is reached at station C. The plateaus be- 
tween stations C and D and between stations E and F represent failure loads of the 
second- and third- stage honeycomb cartridges. If the primary strut begins to extend 
after undergoing compression stroking, the honeycomb load in the strut will reduce 
linearly. This reduction in load is shown in figure 8 along the path from station F to 
station G. If the primary strut continues to extend, the honeycomb load will go to zero 
at station G and will remain at zero if further elongation occurs. If the primary strut 
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starts to compress again, the strut will compress under zero load imtil station G is 
reached, whereupon the load will increase linearly from station G to station F. The 
unloading of the honeycomb from station F to station G is associated with elastic re- 
bo\and of the honeycomb cartridge and may occur at any point in the load-stroke 
diagram. 


Secondary strut 
crushing force, 
FC31 or FC41 


Tension / 




Compression 


The secondary struts are designed such that compression and tension honeycomb 
forces can be generated. A typical secondary strut load-stroke diagram is shown in 
figure 9. The secondary strut load-stroke characteristics are similar to those of the 

primary strut. The undeformed second- 
ary strut starts with zero honeycomb load 
at station zero and strokes in tension or 
compression with a linear buildup in load 
to station A or F. The plateaus defined 
by stations A and B, stations C and D, 
stations F and G, and stations H and I 
represent the failure load of the particular 
honeycomb cartridge. The elastic rebound 
is denoted by the paths from stations D 
to E and from stations I to J. A zero load 
between stations E and J is produced after 
the honeycomb material between stations E 
and J has been crushed. 






stroke. S31 


J or S41 


Figure 9. - Secondary strut load-stroke 
diagram. 


Modification of the Strut Load-Stroke Diagrams 
Because of Structural Elasticity 

The position of the primary and the secondary strut attachment points (or trunnion 
points) on the spacecraft structure are assumed to be fixed in the body coordinate sys- 
tem. In reality, landing gear strut loads are transmitted through these attachment 
points, deforming the spacecraft structure and subsequently causing changes in the posi- 
tion of the attachment points in the body coordinate system. This deformation is im- 
portant in that it effectively softens the elastic character of the honeycomb shock 
absorbers and has the capability to store energy. The influence of the structural de- 
formation on the spacecraft touchdown dynamics can be approximated by idealizing the 
strut- structure system to be as shown in figure 10. This idealization allows for the 

modification of the honeycomb cartridge 
load- stroke diagram to account for the 
structural deformation. 

In figure 10, the spring constant Kg 

is representative of the slope of the load- 
deflection curve of the spacecraft structure 
at the trimnion point and is assumed to be 
constant. The spring constant Kjj is 

representative of the elastic property of 
the honeycomb shock-absorber material. 


K. 


Honeycomb 

5 material 

I— VNA — v^A — m-, 


jn_ 


FVP, FBR 


Figure 10. - Idealized landing gear strut. 
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The terms FVP and FBR are the viscous and friction forces, respectively, present 
during strut stroking. Points A and B in figure 10 are analogous to the end points of 
the struts. 

For the secondary struts, the structural spring constant Kg is made up of two 

factors, the spacecraft structural elasticity and the primary strut bending elasticity. 
The primary stnit bending characteristics are functions of strut length and strut loads, 
which vary during a landing. The primary strut bending and its influence on strut load- 
stroke properties are discussed later in this section. 

The honeycomb cartridge load-stroke diagram is first modified to account for the 
spacecraft structural deformation by computing an equivalent spring rate for springs 
(Ko and K„) in series. The following equation is used to compute the modified onset 

b Xl 

honeycomb spring rates. 

^ ^pS,o 

k^tsKsto 

^ST,o-K^g + Kg,^^^ 

^S^SC,o 

SC,o-Kj,g + Kg^^^ 

modified onset compression spring rate for the primary struts 

modified onset tension spring rate for the secondary struts 

modified onset compression spring rate for the secondary struts 

trunnion-point spring rate at the primary strut- structure attachment 
point 

trunnion-point spring rate at the secondary strut- structure attachment 
point 

primary strut onset spring rate without trunnion-point elasticity 

secondary strut onset spring rate in tension stroking without trunnion- 
point elasticity 

secondary strut onset spring rate in compression stroking without 
trimnion-point elasticity 
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Similarly, the modified primary and secondary strut backoff spring rates are given by 
the following equations. 


_ ^P^P,B ' 

P'B'K^p + Kp B 

_ ^TS^ST, B I 

ST,B=K^S^Kg^ 3 

^ K^rs^sc, B 

sc, B K,j,g + Kg^^ 3 ^ 


(36) 


where KL o = modified backoff compression spring rate for the primary struts 

i', 15 

KgT B “ modified backoff tension spring rate for the secondary struts 

Kgc B ~ modified backoff compression spring rate for the secondary struts 

K_ _ = primary strut backoff spring rate without trunnion -point elasticity 
Jr , r> 

Kst b “ secondary strut backoff spring rate in tension stroking without 
’ trunnion-point elasticity 

Ksc B ~ secondary strut backoff spring rate in compression stroking without 
’ trunnion-point elasticity 

The primary and secondary strut load-stroke diagrams which have been modified 
to include the structural deformation are shown in figure 11. 

The delta displacements shown on figure 11 represent the deformation of the 
spacecraft structure at the trunnion points caused by loading of the primary and second- 
ary struts. The delta displacement of the modified primary strut load-stroke diagram 
is calculated by the equation 


6 


P 


FPl 

Ktp 


(37) 


where FPl is the first-stage crushing force of the honeycomb cartridge in the primary 
struts. Similarly, the delta displacements for the tension side of the modified second- 
ary strut load-stroke diagram are calculated by 


6 


T 


FSTl 

Kts 


(38) 
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(a) Modified primary strut load-stroke diagram. 



Figure 11. - Structural elasticity modification for strut load-stroke diagram. 
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and for the compression side are calculated by 


6 


FSCl 

Kts 


( 39 ) 


where FSTl and FSCl are the first- stage crushing forces of the honeycomb car- 
tridges for tension and compression stroking, respectively, of the secondary struts. 

The modified strut load-stroke diagrams as shown in figure 11 are used as input data 
to the computer. 

A second modification in the secondary strut load-stroke diagram is made to 
accovint for the bending influence of the primary strut. This modification is made in- 
ternally and continuously during the computer simulation. The position of the attach- 
ment point of the secondary struts on the primary strut is assumed to be fixed on the 
idealized rigid primary strut axis. However, because of the primary and secondary 
strut loads on the primary strut, the primary strut will \mdergo bending. The bending 
of the primary strut will cause the secondary strut attachment point to be deflected off 
the idealized rigid primary strut axis. K the assumption is made that the attachment 
point deforms elastically and linearly under load and according to the idealization as 
shown in figure 10, the primary strut bending influence on spacecraft touchdown dynam- 
ics can be approximated by a modification of the secondary strut load-stroke diagram. 

The primary strut is treated as a uniform beam whose elastic axis is defined by 
simple beam theory. The governing differential equation that defines the shape of the 
beam is 


El 


d^?7 


M(s) 


(40) 


where s = position coordinate along the length of the primary strut 
T] = deflection of the center line of the beam 
El = constant bending stiffness 
M(s) = bending moment 

The forces on the primary strut are idealized as shown in figure 12. 

The force F2P is an axial force 
n 

produced by the compression stroking of 
the primary strut. The force FN^^ is the 

normal component of the secondary strut 
honeycomb forces acting on the primary 
strut at the secondary -to-primary strut 
attachment point. This attachment point 


F2Pn 


r\ 


-fNn 

s 


V 

1- 

1 


r n 


-F2Pn 


Figure 12. - Primary strut bending mode. 
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corresponds to s = L21^. By writing the bending moment M(s) as a function of 
F2P^, FN^, and s, equation (40) is expanded to 


^2 / L21 \ 

El ^ + (F2P> = (FN„) (1- V5r^ '<»• 0<=<L21„ («) 


and 


El ^ . (F2P„), = (FN„) (L21J for L21„ < s < V2P„ (42) 

Solving equations (41) and (42) for the deflection t ) and evaluating tj at s = L21^ 
yields 



As the axial force F2P^ approaches zero, the first and second terms of equa- 
tion (43) approach plus and minus infinity, respectively. This situation introduces 
numerical error in the calculation of tj when using equation (43). Therefore, when the 
primary strut compression force F2P^ becomes arbitrarily small, it is assumed to be 

zero and equations (41) and (42) reduce to 


and 



FN 

n 



for 0 < s < L21 

n 


(44) 


El ^ = FN 




for L21 < s < V2P„ 
n n 


(45) 
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Solving equations (44) and (45) for the deflection 77 at s = L21^ yields 


n 


n 



( 46 ) 


Hence, if the primary strut compression force is large, equation (43) is used to com- 
pute the primary strut deflection at s = L21^, and if the primary strut compression 

force is small, equation (46) is used to compute the primary strut deflection. 

The deflection of the primary strut caused by axial and normal forces is assumed 
to be in the direction of the normal force. When a static bending deflection is produced 
in the primary strut at its secondary strut attachment point, an apparent change in the 
length of each idealized secondary strut is produced. This change in length 6. of 
th th 

the i secondary strut in the n landing gear assembly is computed by 


»3,n = \(%- ■^n)l 


(47) 


th 

where = primary strut deflection at point s = L21^ on the n primary 

strut (computed from eqs. (43) or (46)) 

^ 

= unit vector directed along the normal force FN^ 

\^1^, vectors directed along the secondary struts 

th th 

6. = apparent change in length of i ^ secondary strut as a result of n ^ 

primary strut bending 

The deflection of the primary strut as a result of bending caused by secondary 
strut forces is assumed to be a linear function of the secondary strut honeycomb force. 
This assumption along with the idealization presented in figure 10 allows for the calcu- 
lation of an equivalent spring rate for the secondary-to-primary-strut attachment points 
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to account tor primary stimt bending. The spring rates KB4^ for each sec- 

ondary strut V31^ and V41^, which account for primary strut bending, are calculated 
by the following equations. 


KB3„ 


KB4^ 


FC31 


n 


3n 

FC41 

n 


4n 


(48) 


where FC3l^ and FC41j^ are secondary strut honeycomb forces. 

The secondary strut bending spring rates KB3^ and KB4^ change during the 

landing because of landing gear strut articulation and because of varying primary and 
secondary strut shock-absorber loads. For this reason, the bending spring rates and 
the secondary strut load-stroke diagrams are continually updated throughout a landing 
simulation by computing from either equation (43) or equation (46) and then evalu- 
ating equations (47) and (48). 


Viscous Forces 

The mathematical model of the landing gear strut force attenuation mechanism has 
velocity-dependent forces. These forces act in a direction opposite to the strut strok- 
ing velocity and are added algebraically to other axial forces occurring in the struts. 
This arrangement gives considerable flexibility in modeling a landing gear strut shock- 
absorption characteristic. The landing gear strut velocity- dependent forces are de- 
scribed by the equations 


FVP = CP(V2P ) P 
n \ n/ 

FV31^ = CS(V31 J ® ( 

y 

FV41^ = CS(V4lJ ® 


(49) 


where 


th 

FVP = stroking-velocity-dependent force in the n " landing gear 
primary strut 



th, 

FV31 , FV41 = stroking-velocity-dependent forces in the n landing gear 
^ ^ secondary struts 

V2P^, V31^, V41^ = stroking velocities of the primary and secondary struts 

CP, CS = coefficients of the strut stroking velocities used to model a 
given force-velocity characteristic 

y ,y = exponents of the strut stroking velocities used to model a 
^ ® given force-velocity characteristic 

Each landing gear strut is allowed to have individual velocity -dependent force charac- 
teristics defined by data input for the various parameters in equation (49). 


RESOLUTION OF FORCES AT THE FOOTPADS 


Unit vectors which yield the instantaneous direction of each strut in the body coor- 
dinate system have previously been given by equations (28) and (32). If these vinit vec- 
tors are multiplied by the scalar magnitudes of the axial forces determined for each 
strut, the force vectors applied to the footpads by the landing gear can be determined. 


Let F2P be the instantaneous axial force in the primary strut inner cylinder of 
th ^ 

the n gear assembly. Its magnitude is composed of 


F2P = FCP + FVP + FBR 
n n n n 


(50) 


where FCP = force determined from primary strut shock-absorber load-stroke 
characteristics 

FVP^ - primary strut viscous or stoking- velocity -dependent force 

FBR^ = primary strut bearing friction force (developed in the following section) 

Let the instantaneous axial forces in the inner cylinder of the secondary struts be de- 

3k 2k 

noted by F31 and F41 . The directions of these forces are colinear with the previ- 

ously defined vectors V31 and V41 shown in figure 5(b), and their magnitudes are 
given by 


F31 = FC31 + FV31 \ 
n n n 

) (51) 

F41 = FC41 + FV41 
n n n / 
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where FC31^, FC41^ = forces from secondary strut load-stroke characteristics 

FV31^, FV41^ = secondary strut viscous or stroking-velocity-dependent forces 

These forces are applied to the lower end of the primary strut outer cylinder. Fig- 
ure 13(a) shows the axial strut forces. The primary strut bearing friction forces are 
determined from the normal loads on the primary strut bearing surfaces. The method 
for computing bearing loads is developed later in this section. 

Let 


FS = 
n 


F31 +F41 
n n 


(52) 


be the vector sum of the secondary strut forces acting on the primary strut (fig. 13). 
The inertia properties of the landing gear assembly are represented by a mass concen- 
trated at the ball joint of the landing gear footpad. This concentrated mass is referred 
to as the effective footpad mass. The magnitude of the effective footpad mass is deter- 
mined to approximate the inertia characteristics of the landing gear assembly. The 
assumptions involved in this idealization of the mass properties of the landing gear as- 
sembly are discussed in reference 7. 



(a) Axial strut 
forces. 




(b) Forces and mo- (c) 

ments on the 
primary strut 
caused by sec- 
ondary strut 
forces. 


Resolution of 
forces at 
footpad. 


Figure 13. - Forces applied to the footpad. 


The idealization considers the gear struts to be massless so that moments on the 
primary strut are computed without regard to rotational accelerations. Thus 

V2P X FM = X FS (53) 

n n n n 
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Let 


V21 X FS = MS 
n n n 


( 54 ) 


where MS^ is a moment vector perpendicular to the primary strut and to FS^ as 
shown in figure 13(b). Since FM^ is normal to the primary strut at the footpad, the 
modulus of the left-hand side of equation (53) is 


V2P X FM 
n n 


= (V2Pj(FM„)sinf = (V2PjfM^) 


(55) 


and FM^ is given by 


FM 


MS 

n = V^ 


(56) 


The components of FM^ in the body coordinate system can be obtained by constructing 
a imit vector normal to both the primary strut and the moment vector MS^. Thus 



V2P X MS 
'('^2P„)(Ms") 


and 


FM 

n 



FM 

n 


(57) 


(58) 


Bearing Friction Analysis 

The outer cylinder of the primary strut houses the stroking piston and bearings 
as shown in figure 14. The forces Fl^ and F2^ are the normal forces generated at 
the bearings, and their magnitudes depend on the magnitude of FM^. Since Fl^^, F2^, 
and FM^ are coplanar, scalar moment equations for static equilibrium of the inner 
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Figure 14. - Primary strut internal 
bearings. 


cylinder may be written. Summing mom- 
ents about bearing 1 yields 


FM„(m„) - F2„(D2„) = 0 (59) 


where D1 and D2 are the distances 
n n 

defined in figure 14. The length of the 
inner cylinder Dl^^ is a constant, and 

D2^ is determined from an initial bearing 


separation and from the compression stroke S2P^ of the primary strut. The mag- 
nitude of the bearing force F2^ is computed from 




D2 


(60) 


Summing moments about bearing 2 yields 

FM fDl - D2 ) - FI fe2 \ = 0 
n\ n n/ n\ nj 

and 

FM„(D1^ - D2J 

D2„ 


(61) 


(62) 


The magnitude of the bearing friction force FBR^^ can now be calculated from 

(63) 

where i-l„ is the bearing friction coefficient for both bearings. The direction of the 
13 

bearing friction force is taken to be opposite to the stroking direction of the strut. 
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Total Forces and Moments From Landing Gear 

With the bearing friction force now determined, all the terms on the right-hand 
side of equation (53) are known, and subsequently, F2P^ can be obtained from 


F2P 

n 



(64) 


where the unit vector V2P^ is given by equation (28). The total instantaneous force 

transmitted to the footpads by primary and secondary struts (fig. 13(c)) can therefore 
be represented by 


FP = F2P + FM 
n n n 


(65) 


The resultant torque TCMG on the rigid-vehicle center of mass due to landing gear 
strut loads is calculated by the following equation. (For the LM landing system, m = 4. ) 


^ ^ ^ 

TCMG = (bPj^ X -FPj (66) 

n=l' ' 

where BP is the vector from the rigid-vehicle center of mass to the footpad of the 
n ” gear assembly, and -FP is the reaction force applied to the vehicle as a result 
of FPjj- The vector TCMG is expressed in the body coordinate system. 

The translational equations of motion for the rigid-vehicle center of mass are 
written with respect to the inertial coordinate system; therefore, it is necessary to 

express the resultant force fn the inertial system to obtain the inertial acceler- 

ation of the rigid-vehicle center of mass. By using the transformation given by equa- 
tion (3), expressed in the inertial system is 

^ “ - pBl] 
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The sum of the landing gear forces acting on the rigid-vehicle center of mass expressed 
in the inertial coordinate system is then 


m 


FCMG = V FPI 
■^1 n 
n=l 


( 68 ) 


The forces and torques applied to the rigid-vehicle center of mass by the landing 
gear have now been developed. Landing gear forces applied to the footpads have also 
been determined and will be used in the discussion of footpad motion on various landing 
surfaces in a later section. 


ENGINE THRUST AND NOZZLE-CRUSHING FORCES 


Certain spacecraft landing procedures may result in the descent-stage rocket en- 
gine thrust forces being present during the touchdown. Since these forces can be sig- 
nificant, their influence on spacecraft landing performance must be considered. The 
engine thrust forces are usually in a shutdown transient phase at some point in the 
touchdown. The form the thrust tailoff curve will take is a fvinction of the time at which 
the command is given to terminate engine thrust, electrical-mechanical delays in the 
engine systems once the command to terminate engine thrust is given, and the thrust 
tailoff properties of the particular rocket engine. A typical engine thrust tailoff history 
is shown in figure 15. 


The landing procedure used for the 
LM is one where a commanded descent 
rate is given such that the vehicle ap- 
proaches the landing site at a constant ve- 
locity. (In this situation, the rocket thrust 
is equal to the spacecraft weight. ) At 
some preset height above the surface, the 
rocket engine thrust termination command 
is given. This time is represented by tj 

in figure 1 5. The station t^ on the time 

axis represents the position on the thrust 
curve when the footpad touches the surface. 

The thrust curve in figure 15 is approximated in the range from t to t„ by a 
constant Fj^ and from t 2 to t by a polynominal or exponential fxmction. The thrust 
force FTHRUST is given by the equation 


FTHRUST = Fj (69) 



Figure 15. - Rocket engine thrust tailoff 
characteristics. 
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and 


FTHRUST = Fj6 ^ 2 ; ^ ^ ^ + . . . t > tg (70) 

where the constants F^, K, M^, M^, et cetera are determined by fitting the appro- 
priate function to empirically developed rocket engine thrust tailoff data. 

If the location of the rocket engine on the spacecraft is such that the nozzle exit 
plane of the engine comes close to the landing surface while the engine is thrusting, an 
amplification in the engine thrust will occur. The thrust amplification is a complex 
function of such factors as surface geometry, rocket engine design, and engine -surface 
orientation. The amplification factor can be determined by experimental or, in some 
special cases, by analytical means. A typical amplification factor AMP variation with 
nozzle clearance for the LM descent-stage rocket engine thrusting near a flat surface 
is shown in figure 16. The plateau in the thrust amplification curve around zero nozzle 

clearance occurs when the rocket engine 
exit area is less than the throat area. 

The total thrust on the vehicle is 
computed by establishing the distance of 
the centroid of the rocket engine nozzle 
exit plane above the surface (ZNOZ) and 
by determining the corresponding ampli- 
fication factor by evaluating a function 
which describes the characteristics in fig- 
ure 16. The amplification factor is then 
used to multiply the thrust force 

FTHRUST of equations (69) and (70) to 

obtain the total thrust TAF, 



Nozzle clearance. ZNOZ, ft 


Figure 16. - Engine thrust amplification 
factor for the LM. 


TAF = AMP (FTHRUST) 


(71) 




The thrust force TAF is directed through the vehicle center of mass. 

An additional force is involved if the nozzle of the descent-stage rocket engine 
contacts the landing surface. This force is due to the compressive failure forces of 
the rocket nozzle structure. The rocket nozzle structure is assumed to compress 
under constant load as the nozzle structure is forced onto the surface. 

The forces produced on the vehicle because of nozzle crushing are distributed 
around the nozzle exit perimeter in proportion to the nozzle perimeter in contact with 
the landing surface. These forces are directed up the center line of the vehicle. 
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The total force vector due to the engine thrust and nozzle crushing load is 


FpF = TAF + FFL 


(72) 


where FFL is the nozzle crushing load vector. These vectors are referenced to the 



body coordinate system. The moment vector MTF acting on the idealized vehicle be- 
cause of the rocket thrust and nozzle crushing load is 

MTF = Rp X TAF + Rp X FFL (73) 

—V 

where the vector Rp is drawn from the origin of the body coordinate system to the 
centroid of the nozzle exit plane and where the vector Rp is drawn from the origin of 
the body coordinate to the centroid of the nozzle crushing load. The translational force 

S. 

vector FT FI (due to engine thrust and the nozzle crushing load) on the idealized vehi- 
cle in the inertial coordinate system is given by 

FT FI = pgj Fpp (74) 

FORCES AND MOMENTS FROM REACTION CONTROL SYSTEMS 


The forces and moments applied to a landing vehicle by the vehicle RCS are fvinc- 
tions of the particular control system design. The equations used in this analysis for 
determining when control thrusters are activated are simplified approximations to the 
control logic used for the LM. Because of the large attitudes and attitude rates experi- 
enced by the landing vehicle during touchdown, the typical control thruster operation is 
for the thruster to be turned on and to remain on for a large portion of the time interval 
of interest. For this reason, the finer points of control system operation involving 
limit cycles, minimum jet on-time, and onboard computer sampling rates are not of 
prime importance in determining forces and moments applied to the vehicle. Conse- 
quently, relatively simplified control laws are used for landing dynamics simulation. 

Three basic modes of control system operation are considered, including (1) atti- 
tude hold, (2) rate command, and (3) downward translation. The first two modes are 
basic modes used in the LM control system, and the third is a special landing mode 
which has been studied but has not been implemented for the LM. 

It is assumed that the control system applies a positive or negative torque of 
fixed magnitude to the Xg-, Yg-, or Zg-axes, depending on the value of a linear 
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combination of attitude and attitude rate errors associated with the respective body 
axis. The attitude rate errors are computed from 


Eco =0) 

X X X 

Ecd = CO " Cco > 

y y Y { 


Eco„ = co„ “ Cco„ , 
z z z / 


(75) 


where Cco , Cco , and Cco are specified commanded attitude rates. The determina- 
X y z 

tion of attitude errors used in the simplified control laws involves treating attitude 
errors as components of a vector, so that deviations in the values of the Euler angles 
from the commanded values may be transformed into attitude errors which can be com- 
pared with body rate errors. 

The commanded values for the Euler angles are specified by data input and are 
denoted by C0^, and C6^. The deviations in Euler angles from the commanded 

values are computed from 


D0 = 0 

- C0 \ 

X X 

X 

o 

0:5 

II 

o:> 

- C0 > 

y y 

y ' 

D0 = 0 

- ce 


z z z / 


(76) 


The deviations are assumed to be small. With this assumption, these deviations can be 
treated like components of a vector, so that they may be transformed by means of the 
transformation of equation (6) to corresponding attitude errors associated with the body 
system axes. The attitude errors for the body system axes are given by 


E0 " 


D0 1 
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E0 
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D0„ 
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(77) 
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For the attitude-hold mode, the following error equations are evaluated. 


EX . H. ' 

Ey = k2,EeyH.k2jE<„y I 

EZ = kjjES^ + kjjE^^ 


( 78 ) 


where k 


11 ’ 12 ’ 21 ’ 22 ’ 


kgi, and k 32 


are constants, given by data input, which 


model a given control mode and where the commanded attitude rates are zero. The 
errors EX, EY, and EZ are compared to upper and lower bounds given by data input 
to determine if RCS thrusters should fire to produce a torque on the given axis. Let 
EXM, EYM, and EZM be positive numbers given by data input. The torque applied 
to the Xg-axis by RCS thrusters is determined from the following logic. 


-TORQX if EX > EXM ' 
TRCSX=/0 if -EXM < EX < EXM ^ 
+TORQX if EX < -EXM 


(79) 


where TORQX is the magnitude of the torque which can be generated on the X_-axis 

D 

by the RCS thrusters. Similar logic is used to determine TRCSY and TRCSZ, the 
RCS torques about the Y^.- and Z^^-axes, respectively. The RCS torque vector is de- 
fined as 


TRCSX 


MRCS = 


TRCSY 


TRCSZ 


(80) 


The net force vector FRCS due to the RCS would be the null vector for this case. 

The rate-command control mode is a special case of the attitude-hold mode where 
the attitude error bias is removed by setting k^j, kgj, and k^^ equal to zero in 

equations (78) and where the commanded attitude rates are not necessarily zero. Ap- 
propriate bounds EXM, EYM, and EZM for attitude rate errors are specified, and 
the control logic of equation (79) is used to determine control torques. 
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The third control mode, downward translation, is more effective in preventing 
vehicle toppling during touchdown than either of the two modes discussed previously. 
The control logic consists of continuously firing those RCS thrusters which produce 
forces directed downward along the center line of the vehicle. In general, the vehicle 
would be flown in either of the two previously discussed modes for some specified time 
(approximately 1 second for lunar environment) after first contact with the surface; 
then, control would be switched to the downward translation mode. The delay in firing 
downward translation thrusters would avoid contributing to the effective weight of the 
vehicle during the period of maximum energy -absorption requirements on the landing 
gear. 


The downward translation mode is simulated in the landing analysis by specifying 
a delay time from the instant of surface contact for initiation of the downward transla- 
tion. During this delay period, control torques on the vehicle are computed according 
to either of the two previously discussed logics. After the delay period, no pure couple 
torques are applied to the vehicle by the control system, and a force of magnitude 
FRCSZ is applied vertically down the center line of the vehicle. The magnitude of 

FRCSZ is determined from the number of control jets available and from the thrust 



magnitude of each jet. The net force vector FRCS in the body coordinate system as 
a result of this landing mode is 


0 


FRCS = 


-FRCSZ 


( 81 ) 


and then by transformation 


FRCsl 


*Tgjl FRCS 


( 82 ) 


where the vector 


FRCSI 


is the control jet force in the inertial coordinate system. 
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NET FORCE AND TORQUE VECTORS 


The total force vector F (excluding gravity forces) acting on the idealized rigid- 
vehicle center of mass is the sum of individual external force vectors and is given by 
the equation 


F = 


FX 

fy\ 

FZ 


= FCMG + FTFI + FRCSI 


(83) 


Similarly, the net torque on the idealized rigid-vehicle center of mass is given by the 
equation 



T = 


T 

y 


= TCMG + MTF + MRCS 


T 

z 


(84) 


The net vehicle forces FX, FY, and FZ and torques T^, T^, and T^ are used in 

conjunction with equations (10) to (12) and (13) to (15) to calculate the vehicle transla- 
tional and rotational accelerations. 


FORCES ON LANDING GEAR FOOTPADS 


Description of Surface Models 

It has been shown in the preceding sections that, given the position of each land- 
ing gear footpad in the body coordinate system, the instantaneous landing gear geometry 
and landing gear stroking forces can be determined. The position of each landing gear 
footpad is obtained by numerically integrating the equations of motion for each footpad. 

In the following paragraphs, the forces applied to the footpads by the landing surface 
are analyzed, and the equations of motion for the footpads are derived. 

The mathematical surface models developed in this analysis are divided into two 
categories: (1) rigid surface models and (2) nonrigid surface models or soil models. 

A footpad that contacts a rigid surface model is capable only of sliding on the surface, 
and there is no footpad penetration into the surface. However, when a footpad contacts 
a nonrigid surface, the footpad penetrates the surface because of the plastic deformation 
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which the model is capable of simulating. For a given landing simulation, either type 
of landing surface can be specified for any footpad. Consequently, a wide variety 
of surface conditions can be modeled for a touchdown simulation. 


For both surface models, the landing surface coincides with the X-Y plane of the 
inertial coordinate system. This plane is oriented with respect to the gravity vector by 
the transformation of equation (2), which uses an input slope value. Various 

holes and obstacles on the landing surface plane can be specified for each footpad in or- 
der to simulate a given surface topography. .Figure 17 defines the surface irregularity 
parameters HOLE and OBST which are specified for each footpad by data input. 

” The variable HOLE denotes the surface 

n 

elevation of each footpad at the point it 



first strikes the surface. 

For the rigid surface model, sever- 
al types of constraints for lateral motion 
of a footpad along the surface may be 
specified. A positive value for the input 
parameter OBST^ in figure 17 denotes a 

complete lateral constraint on footpad mo- 
tion. This constraint might be visualized 
as a vertical cylindrical wall with its 
height equal to OBST^. Thus, after the 


Figure 17. - Landing surface hole-obstacle footpad contacts the surface, it cannot 

definition. move laterally until it rises above the wall. 

With this option, it is assumed that the 
footpad is rigid and that the sides of the 
wall produce a force equal and opposite to the strut forces resolved at the footpad. K 
a negative value of OBST^ is selected, the footpad is allowed to slide on the surface 

under the influence of coulomb friction, with the friction coefficient specified by the 

variable U . 

n 


Another type of constraint for sliding footpads are curbs which can be placed ar- 
bitrarily on the rigid surface. Two curbs may be specified for each footpad. The curb 
option is specified by using a positive value for the input parameter CLUEW^^. Normal 

forces between the pads and the curbs are developed by using a load-deflection table for 
the footpad. Forces parallel to the curbs are developed by a coefficient of friction UW^ 

and the normal force between the footpad and the curb. The direction of this frictional 
force is considered to be opposite to the sliding velocity of the footpad along the curb. 
Use of the curb option allows for simulation of crushable footpads, since a load- 
deflection table for footpad crushing normal to each curb is specified. When simulating 
landings where the footpads slide along the surface, the bottoms of the footpads are as- 
sumed to be rigid. 


The nonrigid surface model is used to simulate landings on surfaces of various 
bearing strengths. Obstacle and curb options are omitted when this model is selected. 
The simulation of nonrigid models involves equations which are used to generate vertical 
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and lateral forces on a footpad as it translates into the soil. Various forms of these 
equations can be used, and they primarily reflect the variation of bearing strength with 
the position and velocity of a footpad. A detailed discussion of the various landing sur- 
face models and derivations of the footpad equations of motion are given in the following 
sections. 


Rigid Surface Analysis 

In the rigid surface model, there is no movement of a footpad in a direction nor- 
mal to the surface when it is in contact with the surface. Therefore, there is no need 
for an equation of motion for footpad motion in the Z-direction. Likewise, if a foQtpad 
is completely constrained in lateral movement, no equations of motion are needed in 
the X- and Y- directions. 

Since there is no equation of motion to be solved for vertical footpad motion, the 
vertical footpad position is determined through logic in the computer program. At each 
integration step, a tentative Z -position coordinate for each footpad is computed by as- 
suming that the footpad is rigidly connected to the rigid-vehicle center of mass and by 
using landing gear geometry computed in the previous integration step. K this compu- 
tation indicates that the footpad is beneath the local surface plane, the footpad is in sur- 
face contact and the Z-position coordinate is adjusted so that the footpad is on the local 
surface plane. If the computation reveals that the footpad is above the surface and if no 
previous surface contact has been made, then all footpad position coordinates and ve- 
locities are computed as though the footpad were rigidly connected to the vehicle. If 
the footpad is above the surface and if previous surface contact has been made, the 
Z-position coordinate is adjusted so that the footpad is back on the surface as long as 
compressive forces remain in the primary s^trut. (Placement of the footpad back on the 
surface effectively elongates the primary strut and thus unloads the honeycomb shock 
absorbers. ) When all compressive forces in the primary strut have been relieved, the 
tentative Z-position coordinate of the footpad is used as the true footpad position. On 
simulations with lateral footpad constraints, the footpad is not completely free until the 
computed Z-position coordinate of the footpad indicates that the footpad is above the 
constraining obstacle. When the obstacles are cleared, footpad position and velocity 
are computed in the same manner as before surface contact by using newly deformed 
landing gear geometry. 

As mentioned previously, the surface forces on a sliding footpad are determined 
from a coefficient of sliding friction and from forces generated between the footpad and 
a curb which it may contact. It is assumed that frictional forces on the bottom of a 
footpad occur in a direction opposite to the instantaneous sliding direction in the iner- 
tial X-Y plane, and this frictional drag force is computed by the equation 


FDRAGP 

n 



( 85 ) 
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— ^ th 

where FZI = inertial Z-component of the vector FPI on the n footpad (computed 
” in eq. (67)) " 

* ^ th 

SP^ = velocity vector of the n*"" footpad 

th 

= coefficient of sliding friction between the n*^^ footpad and the landing 
surface 

To determine the forces generated between a footpad and a curb, it is necessary 
to describe the lateral crushing characteristics of the footpad by data input. A typical 
crushing characteristic for a footpad is shown in figure 18, 

Curbs are oriented on the surface with respect to the initial footpad contact point. 
Two curbs can be designated for each footpad by the parameters shown in figure 19. 



Crush distance normal to curb, ft 

Figure 18. - Force-deflection diagram 
for footpad crushing. 


The inertial coordinates of a footpad 
are used to compute the distance from the 
center of the footpad to a curb at any time 
while the footpad is on the surface. If this 
distance is less than the original footpad 
radius, the amount of footpad crushing 
normal to the curb can be determined. 

The force between a footpad and a curb is 
then obtained from the footpad force- 
deflection diagram in figure 18. 

When a footpad first contacts 
the surface, the inertial position coor- 
dinates of the footpad are denoted by 

rxPI , YPI , -HOLE These coordi- 
n n n-^ 

nates are determined by transforming 
footpad coordinates in the body system 
to inertial coordinates. The coordinates 



DWALLI . = the normal distance from 

^ th th 

n footpad center to j 

curb when footpad initially 

contacts the surface 

WALANG . = the angle from the positive 
^ X-axis to a line in the 
X-Y plane normal to the 

curb of the n^^ footpad 

Figure 19. - Schematic diagram of curb 
parameters on rigid surface. 
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of position vectors W which locate the intersection point of DWALLI with a curb 

n, j n, j 

are computed from 


W . = f/XPI +a A,/YPI +b 

n, j ” n, j; ’ n n, ’ J 


( 86 ) 


where 


a . = DWALLI . cos /WALANG 
n, J n, J V n, j) 


b . = DWALLI . sin /WALANG A 
n, J n, ] ( n, j) 


(87) 

( 88 ) 


Figure 20 illustrates a typical initial footpad-curb position. 



Figure 20. - Initial footpad-curb position 
on rigid surface. 


Then, a unit vector 6 . lying along the 

o, J 

curb can be constructed by the vector 
cross product 



= . x[0, O,!]”^ 


(90) 
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At any time after a footpad initially contacts the surface, vector A . joining the foot- 

V J 

pad center with the point located by W . can be constiructed so that 

n, 1 


A . = 
n,3 


XPI + a„ . - XP. 
n n, ] 


n 


YPI + b„ . - YP„ 
n n, ] n 


(91) 


where XP and YP are the instantaneous inertial coordinates for the footpad center, 
n n ^ 

The current distance between the footpad center and the curb can then be computed from 

the vector dot product 


DWALL . = A . • AI . 
n, 3 n, 3 n, 3 




DWALLIn^.^, 


?XPI + a . - XP Va .) 
I, n n, 3 nj\ n, jJ 


(92) 






by 


Given the undeformed footpad radius PR, the variable PSTRK^ ^ is computed 


PSTRK . = PR - DWALL . 

j ^>3 


(93) 


If PSTRK . is a negative quantity, the footpad is not in contact with a curb; but if 
3 

PSTRK . is positive, its value represents the amount of footpad crushing normal to 
th ^ 

the 3 curb. Appropriate logic is used to retain deformed footpad geometry for multi- 
ple impacts between footpads and curbs. 


A crushing force value PFRC . is obtained from the footpad force-deflection 

3 

diagram (fig. 18) by using PSTRK .. Since PFRC . acts on the footpad in a dir ec- 

3 3 

tion normal to the curb, PFRC . can be written as the force vector 

’ n,3 


PFRC„ . 
n, 3 




AI 


n,3 


(94) 
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other forces on the footpads are generated at the curbs by friction between the 
footpad and curb. These frictional forces are computed from PFRC^ ^ and a coeffi- 
cient of friction UW^. Frictional forces are applied to the footpad parallel to the curb 
but opposite to the velocity of the footpad along the curb. A scalar product of SP and 
the unit vector ^ along the curb is performed to obtain the proper sign for the fric- 

on a footpad is then given by 


tion force. The curb frictional force FWF 


n,] 


FWF . = ±(^UWYpFRC .'iB . 
n,j \ n,]; n, ] 


(95) 


It is assumed that the direction of 
gration step. 


FWF . 
n,] 


does not change during a numerical inte- 


All the instantaneous forces FPI^, FORAG^^, PFRC^^, and FWF^ acting on a 

footpad are now accounted for, and the instantaneous vector equation of motion of a 
footpad is obtained from a vector sum of all forces acting on the footpad. 


SP = 
n 


^FPI + PFRC + FWF ^ 
n n n 

, PMASS y 


y^FDRAGP \ SPj^ 


PMASS 


SP 


+ G 


(96) 


where PMASS = effective mass of a footpad 

G = gravity vector in inertial coordinates 

For computation purposes, equation (96) must be reduced to component form. Then, 
the equation becomes 


FXP 


XP = 


n 


FDRAGP XP 
n n 


n PMASS 


PMASS /(XPJ2 ^ (y-p J 


FYP 


YP = 


n 


FDRAGP YP 
n n 


n PMASS 


PMASS /(XPJ -.(YP^) 


2 y 


/ 


(97) 


ZP = 0 
n 
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where 


FXPn 

FYP„ 


r FPI + PFRC + FWF 1 
L n n nj 


Fm + PFRC + FWF 1 
_ n n nj 


[1,0,0]'^ 

(98) 

[0,1,0]'^ 

(99) 


and = component of gravity acceleration in the inertial system X- direction 
gy = component of gravity acceleration in the inertial system Y-direction 

Equations (97) are integrated by a Rvinge-Kutta method, where FXP^, FYP^, 

and FDRAGPj^ have constant values over the integration interval. However, because 

of changes in footpad velocity components during the integration step, this method of 
solution does account for changes in direction of the drag force. Equations (97) have a 

singularity if the footpad velocity SP^^ becomes zero. Since the equations are solved 

numerically over discrete time intervals, it is unlikely that the net result of any inte- 
gration time step would yield a velocity value of exactly zero at the end of the interval. 
Instead, the numerical integration would actually change the direction of the velocity 
component. Such an occurrence would imply that the frictional force is capable of de- 
celerating the footpad to zero velocity and then accelerating it in the opposite direction 
during the integration interval. Since in reality, a frictional force should not be able 
to accelerate a motionless footpad, the following special logic is used to avoid this 
problem. 

If the sign of either XP^^ or YP^^ is reversed during any integration step, it is 

assumed that the particular component of velocity is integrated to zero during the in- 
terval. The position of the footpad is then adjusted by determining an average velocity 

for the component during the integration interval. If the signs of both XP^ and YP^ 

are changed during an integration time interval, then the footpad has stopped, and it 
is assumed that it remains at rest until the forces applied to it by the landing gear and 
the curbs are sufficient to overcome the friction force. The vector which represents 
the sum of the X- and Y-forces applied to the footpad by the landing gear and curbs is 
defined by 

FXYPj^ = [ FXP, FYP, 0] (100) 


The magnitude of FXYP^^ is compared with FDRAGPj^. Whenever FXYP^^ is larger 
than FDRAGP^, one integration step is made, with X- and Y-components of FDRAGPj^ 
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determined by assuming FDRAGP^ acts in a direction opposite to FXYP^^. After this 

integration step, nonzero values for XP and YP are obtained, and equations (97) 
can again be used. 

Numerical integration of the footpad equations of motion yields values for each 
footpad coordinate in the inertial system. These coordinates are transformed back into 
the body system, thus allowing for computation of new landing gear strut geometry and 
landing gear forces for the next integration interval. 


Nonrigid Surface Analysis 

The physics of the response of soil media to projectile impacts is not a well- 
developed science. Until recently, the best available knowledge of forces produced in 
soils could be obtained only from texts oriented to the study of building foundations 
(ref. 9). An interest in the ability of spacecraft to land on the lunar surface and spec- 
ulation about the bearing properties of the lunar surface have fostered a considerable 
amovmt of research devoted to understanding the soil impact phenomena (refs. 6 and 
10 to 12). Most of this research has been experimental, and only empirical methods 
of predicting forces which are applied to an impacting body by a soil medium have been 
developed. These methods generally attempt to relate soil forces on the impacting body 
with the penetration velocity of the body, along with the physical properties of the soil 
(such as density, relative density, cohesion, internal angle of friction, etc. ). To gain 
some insight into the dynamic behavior of vehicles landing in various types of soils, 
mathematical soil models are useful. 

Although several different mathematical soil models have been used in spacecraft 
landing simulations, only one of these models is discussed in this report. In this 
model, no attempt is made to relate the forces generated by the soil to the physical 
properties of the soil. Instead, the soil forces are related to footpad penetration and 
footpad velocity by constants, which represent penetration -resistance characteristics 
of the soil. 

Forces in the inertial Z-direction, exerted on the footpad by the soil, are com- 
puted from an equation of the form 


FZP = AP 
n 


K,*K2(zpEnJ.K3(zpJ.K,(zpJ2 


where AP = area of the bottom of a footpad 

Kj = bearing strength of the soil at the surface 
K 2 = increase in soil bearing strength with depth 
ZPENj^ = penetration of a footpad in the inertial Z-direction 
Kg = increase in soil bearing strength with footpad velocity 


( 101 ) 
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= increase in soil bearing strength with the square of the footpad velocity 
ZP^ = penetration velocity of a footpad 

The forces applied to the footpad in the X-Y plane are obtained from a similar equation 


FDRAG = -APF 
n 


■^5 ♦ * sKf 


SP 
n 

SP 


( 102 ) 


where APF = frontal area of a footpad in contact with soil (assuming the footpad is 
cylindrically shaped, and computing area as a function of penetration) 

Kg = static soil resistance to horizontal movement 

Kg = increase in soil resistance to horizontal motion with respect to footpad 
lateral displacement 

SPj^ = lateral translation of a footpad 

K„ = increase in lateral resistance of the soil with respect to the depth of 
penetration of the footpad 

Kg = increase in soil resistance with respect to the square of the footpad 
horizontal velocity 

SP^ = instantaneous horizontal velocity of a footpad 

By changing the constants in equations (101) and (102), numerous soil characteristics 
may be simulated. If lunar and planetary landings are to be studied, equations (101) 
and (102) can be as useful as some of the more complex soil impact models 
(ref. 6), since soil properties at a particular landing site will not be well defined and 
since it is desirable to simulate surface strength characteristics of a widely varying 
nature. 

The equation of motion for inertial Z-translation of a footpad is obtained from the 
summation of forces in the Z -direction 


ZP 


= /'fZP - FZI 


n V 


n 


n) 


PMASS ^z 


(103) 


Equation (103) is solved numerically by assuming constant forces over an integration 
interval. The X and Y translation equations of motion are identical to equations (97), 



since FDRAG has the same form as FDRAGP in these equations, 
n n ^ 
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A new variable FGRND is defined by 

n 


FGRND = FDRAG + FDRAGP 
n n n 


(104) 


and is substituted for FDRAGP^^ in equations (97). This substitution allows for fric- 
tional forces on the underside of the footpads to be added to the soil forces of equa- 
tion (102). Equations (97) are integrated in the same manner described previously. 


STABILITY CRITERIA 


A performance evaluation of a soft- landing system requires the simulation of 
many touchdown conditions. When a large number of touchdown simulations are re- 
quired, the speed at which the computer can execute these simulations becomes a fac- 
tor in the problem. In the evaluation of a soft-landing system for stability (tipover) 
performance, the computer running time can be reduced by defining stability criteria 
such that the computer can make the decision as to whether the landing being simulated 
will eventually be stable or unstable. Upon making this decision, the computer simula- 
tion is terminated. 

A parameter representative of the state of stability of a particular landing is the 
magnitude of the minimum "stability distance." The stability distance is the shortest 
distance between the rigid-vehicle center of mass and a plane parallel to the gravity 
vector that passes through two adjacent landing gear footpads. This plane is defined 
as the "stability wall. " Presented in figure 21 is the vector definition of the stability 

distance. The length of vector SD is the 
stability distance. The stability distance 
SD is calculated by the equation 


gp _ ACMj ((^ AB) 
Gx AbI 


where AB = vector connecting two adja- 
cent landing gear footpads 


ACM = vector connecting the rigid- 
vehicle center of mass and 
a landing gear footpad 

G = gravity vector 



Figure 21. - Vector definition of the 
stability distance. 
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The stability distance vector is calculated by the equation 


^ SD(G X AB) 
G X AB 


( 106 ) 


A stability distance is computed in this manner for each stability wall. 

If the rigid-vehicle center of mass is within an enclosure defined by stability 
walls passing through adjacent landing gear footpads, then the gravity action on the ve- 
hicle dynamics is stabilizing. (The center of mass and center of gravity are assumed 
to be coincident. ) If the rigid-vehicle center of mass is outside of this enclosure, the 
gravity action on the vehicle dynamics is destabilizing, and the landing is declared un- 
stable. If the vehicle is unstable, the degree of instability can be gaged by the magni- 
tude of the "overturning velocity" at the instant the rigid-vehicle center of mass passes 
through the stability wall. The overturning velocity V is calculated by the equation 


V = V ■ ^ 
o 


(107) 


where V is the translational velocity vector of the rigid-vehicle center of mass. 

Also, instantaneous values of the stability distance, overturning velocity, and 
vehicle kinetic energy can be used in determining if the landing will eventually be sta- 
ble. At the instant a landing gear footpad contacts the landing surface, the vehicle has 
a given kinetic energy. Subsequent motion will increase or decrease the initial vehicle 
kinetic energy, depending on the amount of energy dissipated by the landing gear sys- 
tem or the energy gain due to the lowering of the rigid-vehicle center of mass. When 
the vehicle kinetic energy has been dissipated to the point that it is less than approxi- 
mately 10 percent of the initial vehicle kinetic energy at touchdown, the landing gear 
struts have usually undergone maximum stroking. After maximum stroking of the 
landing gear struts occurs, the resulting vehicle motion will generally be rigid-body 
rotation about two adjacent landing gear footpads on the landing surface. .After com- 
puting the minimum stability distance to each stability wall, a landing can be declared 
stable if it satisfies the following conditions: 

1. The vehicle kinetic energy is less than 10 percent of the initial touchdown 
kinetic energy. (This constraint can be replaced by not declaring a landing to be stable 
until landing gear struts have completed stroking. ) 

2. The rigid-vehicle center of mass is inside the stability walls. 

3. The two adjacent landing gear footpads for which the minimum stability dis- 
tance is computed are on the landing surface. 

4. The rigid-vehicle center of mass is moving away from the stability wall with 
some preselected minimum velocity. 
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An index of the degree to which the landing was stable can be judged by the magnitude 
of the minimum stability distance obtained during the landing. 

The capability of the computer program to decide whether a landing is stable or 
imstable enables the program to be operated in a mode where a "stability boundary" 
can be defined with a minimum number of landing simulations. The stability boundary 
is a curve of vertical touchdown velocity compared to horizontal touchdown veloc- 

ity Vj^; this curve separates the plane into regions of velocity conditions that 

result in either stable or unstable landings.. The generation of the stability boundary 
is accomplished through program logic, which enables the program to choose touch- 
down velocity conditions to be simulated, based on the stability results of previous 
landing simulations in the same computer run. This capability is of considerable im- 
portance when extensive studies of landing stability performance are required. 


ANALYTICAL-EXPERIMENTAL DATA CORRELATION 


The degree to which the mathematical model of a soft- landing system represents 
physical reality depends primarily upon the assumptions made in the analysis. In eval- 
uating landing performance of a soft-landing system, factors such as vehicle stability 
and landing gear energy-absorption characteristics are of prime importance. In light 
of the complexities of a soft- landing system, some approximations are inevitable, and 
for this reason, experimental verification of the analysis is desirable. 

Several experimental programs have been conducted to develop landing dynamics 
data. Various LM models have been drop tested at the Langley Research Center, at 
MSC, and at NASA contractor facilities. References 13 to 15 present some of the data 
from these test programs. Reference 15 contains comparisons of experimental results 
of full-scale LM drop tests and theoretical predictions of the mathematical analysis 
presented in this report. 

The analytical-experimental data correlation presented in this report is restric- 
ted to landing dynamics simulations of a 1/6-scale LM model. The experimental data 
were developed by a NASA contractor (ref. 13) by dropping a l/6-scale LM model on a 
plywood surface (fig. 22). Critical stability and landing gear energy-absorption data 
were obtained by varying the model touchdown conditions. 

The forces acting on the 1/6-scale model during landing are gravity forces and 
landing gear strut forces. The landing gear strut forces are due to friction and to the 
deformation of the honeycomb cartridges housed inside the struts. The landing gear 
footpads are spikes that penetrate into the plywood landing surface and constrain the 
footpad from lateral motion during impact. The effects of the landing gear footpad 
penetration into the landing surface on the model touchdown dynamics are accounted 
for by modifying the primary strut load-stroke diagram. 

A schematic drawing of the 1/6 -scale LM model is presented in figure 23. 
Landing gear strut load -stroke diagrams are presented in figure 24. 
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Force, lb Force, lb 




(a) Side view. 


(b) Top view. 


Figure 23. - The 1/6-scale LM model geometry. 



(a) Primary strut load- stroke (b) Secondary strut load- stroke 

diagram for a l/6-scale diagram for a 1/6-scale 

LM model. LM model. 



(c) Landing surface load- 
deformation charac- 
teristics. 

Figure 24. - Model physical 
characteristics. 


The geometric parameters and sur- 
face deformation properties for the model 
are as follows: 


R1 = 1.639 ft 
R2 = 1. 562 ft 
R3 = 2. 221 ft 


/3 =19.72° „ 

I = 0.994 slug -ft 

Z Z c 

1 ,I = 1.397 slug -ft" 
xx’ yy ^ « 

I ,I ,I = 0 slug-fr 
xy xz yz 


Z1 = 0.490 ft 
Z2 = 0. 806 ft 
Z3 = 1.351 ft 
ZC = 0.952 ft 


Bearing friction coefficient = 0. 10 
Bearing overlap = 0. 396 ft 
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Two model touchdown orientations were tested and simulated. These touchdown 
orientations resulted in planar motion of the model. One orientation is the planar 2-2 
landing mode. This type of landing occurs when the two trailing landing gear footpads 
strike the surface simultaneously, followed by the simultaneous impact of the two lead- 
ing footpads. 

The second orientation is the planar 1-2-1 landing mode. This type of landing 
occurs when the trailing landing gear footpad contacts the landing surface and is fol- 
lowed either by simultaneous contact of the two midgear footpads or by contact of the 
leading footpad. The 2-2 landings pose critical problems for stability and for second- 
ary strut compression stroking. The 1-2-1 landings present critical primary strut 
compression stroke-out problems. 

Figure 25 shows a comparison of the analytically predicted and the experimen- 
tally measured model vertical acceleration during the first 0. 16 second of a landing. 
The l/6-scale LM model was oriented such that a planar 2-2 landing occurred. After 
the first 0. 12 second of the landing, the resulting motion of the model is rigid-body 
rotation about the two downhill landing gear footpads. The model was initially pitched 
up 5° with respect to the horizontal, with horizontal velocity in the downhill direction. 
The landing surface was sloped 5°, with 4-inch depressions superimposed on the slope 
under the downhill landing gear footpads. 



Figure 25. - Analytical-experimental correlations of the vertical acceleration 

of a 1/6 -scale LM model. 
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The major difference between the analytical and experimental time histories 
shown in figure 25 is the spiked oscillation of the experimental curve. The cause of 
this phenomenon has not been established, but a likely reason may be that excessive 
binding of the primary strut inner cylinder occurred during stroking. Some binding is 
expected and is accmmted for in the analysis by bearing friction forces. However, a 
constant bearing friction coefficient is assumed, and possible deviations from this 
idealization could occur if manufacturing tolerances of the model landing gear are not 
equivalent to prototype hardware. Some wearing of the bearing surface could also 
occur after several drops in the drop-test series. 

The following explanation of the various features of the acceleration time history 
in figure 25 is made with reference to the analytical curve. The first acceleration 
pulse shown in figure 25 is due to the initial impact of the trailing landing gear footpads. 
Subsequent model rotational and translational motion brings the downhill landing gear 
footpads in contact with the surface, which produces the second acceleration pulse. 

After the first acceleration pulse, the model goes into a free-flight period of approxi- 
mately 0. 02 second. The first acceleration pulse is produced predominantly by the 
first-stage honeycomb-crushing forces in the trailing landing gear primary struts. The 
second acceleration pulse is predominantly due to the primary strut loads of the down- 
hill landing gears. The first plateau of the second acceleration pulse is mainly due to 
the first- stage honeycomb- crushing force of the primary strut. The second plateau in 
the second acceleration pulse is due to the secondary strut forces acting on the primary 
strut, producing an increase in the primary strut bearing friction force. The final peak 
in the acceleration pulse is due to the primary strut stroking into the second stage of 
the honeycomb shock absorber, hi this case, the primary struts stroked a distance of 
0. 01 foot into the second-stage honeycomb shock absorber. Shown in figure 26 is a 
comparison of analytical and experimental model pitch rates for the same landing. The 



Figure 26. - Analytical- experimental pitch rate correlation of a 1/6-scale LM model. 
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initial rise in pitch rate is predominantly due to the forces of the trailing landing gear 
primary struts. As the model rotates, the vertical forces are diminished to the point 
that the lateral constraining forces produce a positive pitching acceleration before the 
model goes into the period of free flight. At approximately 0. 06 second, the downhill 
landing gear footpads strike the landing surface. The resulting forces reduce the pitch 
rate from -5. 5 rad/sec to -2. 0 rad/sec. After 0. 12 second, the model has experienced 
maximum stroking in the landing gear struts, with the resulting motion being near- 
rigid-body rotation about the two forward footpads. 


The model drop-test data used in the detailed correlation shown in figures 25 
and 26 were the only detailed time-history data available from this particular model 
test series, and the data were not selected because they represented the best analytical- 
experimental correlation. In fact, these data probably represent one of the poorest 
correlations expected, since the model did not overturn in the test, while the computer- 
simulated landing did overturn. As shown in figure 27, the computer program correctly 

predicted stability results for the other 
model drops in the test series. The 
touchdown attitudes, surface slope, and 
surface depressions for the test results 
in figure 27 are the same as those of the 
drop-test data in figures 25 and 26. 


» 8 
8 


CD 7 


UV\ stability profile: 

2-2 landing 

5“ slope, 5° nose -up 
Bearing friction coefficient = 
4-inch depressions at 


0. 10 


forward gear 



The predicted stability boundary is 
within 1 ft/sec of the experimental bound- 
ary at the vertical touchdown velocity of 
10 ft/sec. For other values of vertical 
velocity, the analytical profile is within 
the bounds of the experimental profile. 

The overall correlation between analyti- 
cal and experimental stability results 
seems good. The data indicate that the 
stability profiles for an LM-type landing 
gear system can be computed to within 
1 ft/sec for landings on rigid surfaces 
with fully constrained landing gear foot- 
pads. Correlations of computer simula- 
tions and drop tests for maximum primary 
and secondary compression strokes are 
presented in table I. 


2 _ O Experimentally stable 
• Experimentally unstable 
I - MSC analytical stability boundary 

I 1 l I 1 i I ! 1 

0 1 2 3 4 5 6 7 8 9 

Horizontal velocity, ft/sec 

Figure 27.- Anal 3 d;ical- experimental 
stability correlation of a 1/6- scale 
LM model. 


The overall correlation between the 
predicted and the measured parameters 
of the l/6-scale LM model landings was 
good because the load- stroke properties 
of the landing gear struts and the landing 
surface force-deformation properties 
were well defined. If the same degree of 
definition can be obtained for a full-scale 
landing vehicle, the predicted landing per- 
formance would have comparable accuracy. 
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TABLE I. - CORRELATION OP ANALYTICAL AND EXPERIMENTAL 


LANDING GEAR STRUT STROKE DATA 


Landing 

Touchdown velocity, 
ft/sec 

Amount of stroke, in. 

(a) 

V 

V 


Experimental data 

1 

Analytical data 



Primary strut 

- 1 ■ 


1-2-1 

-5. 0 

7.0 

3,81 

3.934 

1-2-1 

-7. 0 

6.2 

3.88 

4.264 

1-2-1 

-10.0 

4.0 

4.07 

4.050 

1-2-1 

-13.0 

1.2 

3,91 

3.988 



Secondary strut 


2-2 

-4.0 

7.0 

1. 16 

^ 1. 007 

2-2 

-6.0 

8.0 

1. 86 

1.804 

2-2 

-8.0 

6.0 

2. 33 

2. 358 

2-2 

-10.0 

6.0 

2.72 

2.771 


^1-2-1 landing: 5° nose-up; downhill landing; 5° surface slope; 4-inch depres- 
sions under midgear and downhill gears. 2-2 landing: 5° nose-down; uphill landing; 

5° surface slope; 4-inch depressions under trailing gears. 


CONCLUDING REA/IARKS 


The equations governing the touchdown dynamics of a soft-landing spacecraft have 
been derived and presented. A digital computer simulation based on these equations 
has been developed. The computer program has been utilized successfully in the de- 
velopment of the Apollo lunar module landing gear system and for predicting the lunar 
landing performance of the lunar module. 

The accuracy of the computer simulation of the lunar module landing was estab- 
lished through correlation with data measured during tests of a dynamically scaled 
model. Comparison of data indicates that touchdown velocities which result in space- 
craft overturning can be predicted to within 1 ft/sec. Landing gear strut stroking can 
be predicted to within 10 percent, provided significant stroking occurs. 
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Theoretical methods have proved to be satisfactory for guiding landing gear sys- 
tem development and for determining spacecraft landing performance. This fact, cou- 
pled with the problems associated with providing an adequate test environment, will 
undoubtedly result in more emphasis being placed on theoretical procedures, such as 
the ones presented in this paper, for future spacecraft landing system development. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, September 4, 1970 
914-50-20-89-72 
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APPENDIX 


PROPELLANT SLOSHING 


It Is likely that future spacecraft with a soft- landing capability, like the Apollo 
LM, will use liquid propellants in their primary propulsion systems. During the de- 
scent and subsequent touchdown of the spacecraft, liquid propellant is consumed, re- 
sulting in partially filled propellant tanks at the instant of touchdown. Upon initial 
spacecraft contact with the landing surface, the residue propellant is subjected to a 
sloshing phenomenon which influences the touchdown performance of the spacecraft 
landing system. 

The mathematical model used to simulate the propellant sloshing during the 
spacecraft touchdown is based on ’’linear sloshing theory” (refs. 16 and 17). Linear 
sloshing .theory is used to compute the sloshing force on the side of the propellant tank. 
If the assumption is made that the sloshing force on the tank walls is due to the first 
sloshing mode (as defined in refs. 16 and 17 and fig. 28) and that the higher sloshing 
modes are negligible, then the sloshing forces can be calculated by 


F = F Z„q 
w s 


(108) 



Figure 28. - Sloshing force distribution 
on propellant tank wall (first sloshing 
mode). 


where F = total force applied to the tank 
wall by the sloshing propel- 
lant, lb 

Fg = force per xmit vertical tank 
acceleration caused by unit 
sloshing amplitude, 

(lb/ft)/(ft/sec^) 

q = sloshing amplitude, ft 

z" = vertical acceleration of the 
propellant tank in the body 

2 

coordinate system, ft/sec 


The force F is determined from linear sloshing theory (ref. 18) by integrating the 
s 

force distribution on the propellant tank wall due to a unit sloshing amplitude. Fig- 
ure 28 illustrates a typical sloshing force distribution on the tank wall. The force F 

s 

is located at the centroid of the force distribution, a distance £ below the undisturbed 
free surface level. 
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The equation which governs the sloshing amplitude is given by 


q = 



w 



(109) 


where q„ = rigid-vehicle acceleration of the propellant tank normal to the vertical 

2 

axis of the tank, ft/sec 

2 2 
CO = sloshing frequency squared per unit vertical tank acceleration (rad/ sec) / 

(ft/sec^) 

Zg = vertical acceleration of the propellant tank 

The frequency term co is determined from linear sloshing theory (ref. 18) and is a 
function of such factors as propellant tank shape, propellant loading, and gravity field. 
Numerical integration of equation (109) gives the sloshing velocity and amplitude-time 
histories which are used in equation (108) to calculate F^. 

When a partially filled propellant tank is subjected to a constant lateral accelera- 
tion, the inertial resistance of the sloshing propellant is reflected as a force on the side 
of the propellant tank. Consider the residue propellant mass as being comprised of a 

fixed mass M„ and a sloshing mass M . When residue propellant mass is considered 
1 s 

to be two masses, a linear spring- mass system can be constructed which is analogous 
to the sloshing dynamics. 

In the linear spring- mass system, the force on the propellant tank wall as a re- 
sult of sloshing is 


F = M w q 
w s n ^ 


( 110 ) 


where M = sloshing mass, slugs 
s 

= natural sloshing frequency, rad/sec 

2 2 " 

By equating equations (108) and (110) and by noting that = co then 
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and 


M 




( 112 ) 


Once the sloshing mass is determined, modifications to the total spacecraft mass char- 
acteristics must be made in order to account for the separation of the sloshing mass 
from the remaining vehicle mass. 

Figure 29 illustrates the location of the propellant masses and the rigid-vehicle 
center of mass, with and without propellant sloshing mass. The propellant tank is as- 
sumed to be symmetric about the Z_-axis. 

r> 

The center of mass of the sloshing pro- 
pellant Mg is distance £ below the 

quiescent propellant surface, as is the 
centroid of the sloshing force F^. By 

taking moments about the center of mass 
of the rigid vehicle with sloshing propel- 
lant, II 2 is obtained from 


M 


where M is the total mass of the space- 
craft, including the propellant 
mass. 

At this point, the idealized rigid 
spacecraft is considered to be a two- 
mass system; the two independent masses 
are the mass of the rigid vehicle exclud- 
ing the sloshing mass (M - M ) and the sloshing mass itself. The influence of one 

s 

mass on the other is the propellant sloshing force. Mass characteristics for landing 
vehicles are most often available in a form which assumes the propellant to be a rigid 
part of the spacecraft. To properly incorporate this type of data into landing dynam- 
ics simulations which consider propellant sloshing, the mass data require some 
modification. 

Modification of the vehicle mass moments of inertia involves only those mass 
moments about the Xg- and Yg-axes, because the sloshing mass is allowed essen- 
tially two degrees of freedom (motion in the Xrj-Y„ plane). The change in the mass 

D JD 


ir 

tr- 


I dealized-rig id -vehicle 
(spacecraft) center of mass, 
without sloshing mass 

M, rigid-vehicle (spacecraft) 
center of mass, with 
sloshing mass 



center of mass of 
sloshing propellant 


Quiescent propellant 
surface 


M^, center of mass of 


fixed propellant 


Figure 29. - Influence of sloshing pro- 
pellant on the rigid-vehicle center 
of mass. 
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moment of inertia about the Z„-axis is assumed negligible because of small-amplitude 

Jc3 

displacements. The following equation gives the new mass moments of inertia for a 
rigid spacecraft, corrected for removal of the sloshing mass. 




/ 

\ 2 / 

^ > 

,2 


I 4. 

= I - 

(M - 

M - M fJ 

+ f 



XX, f 

XX 


s) 2 s\^ 

1 > 
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I f 

= I - 

(m - 

M -M (i 

e, + 


> 
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yy 

\ 

sj 2 s\ 
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^zz, f 

= 'zz 





J 


(114) 


where I ,, I I „ = modified mass moments of inertia in the body system of the 
XX, yy, zz, rigid vehicle with the sloshing mass excluded 

I , I , I = mass moments of inertia in the body system of the rigid ve- 

XX yy zz j^icle, including the total residue propellant mass 

The sloshing force vector F^, which is the interacting force between (M - M^) 
and M , has body system components given by 


F = 
w 


PB,x 


PB,y 


(115) 


The body system components are transformed into inertial coordinate system compo- 
nents by 




F„„ 

B, X 


PB, X 

F 

B,y 

II 

1 1 

W 

HH 

1 1 

F 

PB,y 



0 

B, zj 




(116) 


where ~ transformation matrix given by equation (3) 

Ft 5 , F„ , F_ = sloshing force components in the inertial coordinate system 
B, X B, y B, z 
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The sloshing forces create torques about the center of mass of the idealized rigid 

vehicle. The torques T. and T„ about the X„- and Y„-axes are calculated from 

I, X I, y n 


"^f,x ^BP,y (^1 "^^2 I 

"^f,y "" ^BP,x (^1 ^^^2 j 


( 117 ) 


The sloshing forces and torques are considered to be external forcing functions in the 
solutions of the rigid-vehicle equations of motion. 
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